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SCATTERING  OP  A  SPHERICAL  WAVE  ON 


TRUNCATED 


REVOLUTION 


RORIES 


Yu.  A.  Yerukhimovich 


Summary 


A  scattering  field  is  analyzed  on  truncated  revolution  bodies 
formed  by  curves  of  the  second  order  when  they  are  irradiated  from 
one  focus  by  a  source  of  a  spherical  wave.  The  analysis  uses  the 
method  of  curvilinear  coordinates  II).  The  solution  is  obtained  in 
the  distant  zone,  in  approximation  of  physical  optics.  The  system 
of  formulas  describes  the  diagram  of  the  scattering  field  in  all 
space,  including  the  transitional  zones  between  the  regions  of  the 
angles  of  the  observation  point. 
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Int-  duct  ion 

'ublication  {1]  presents  a  method  for  solving  diffraction  problems 
based  on  the  introduction  of  a  special  system  of  curvilinear  coordinates 
associated  with  the  physical  features  of  the  phase  structure  of  the  field. 
In  this  case  ,  the  center  of  the  emitted  spherical  wave  and  the  obser¬ 
vation  point  serve  as  the  focal  points  for  families  of  ellipsoids  and 
hyperboloids  which  form,  two  systems  of  coordinate  surfaces.  The  third 
system  of  coordinate  surfaces  is  half  planes  passing  through  a  line 
which  connects  the  focus  and  the  observation  point.  The  integral 
presentation  of  the  diagram  of  the  scattering  field  in  this  case 
acquires  the  simplest  appearance. 

In  order  to  obtain  a  strict  solution,  it  is  necessary  to  know  the 
exact  quantity  of  the  current  which  flows  on  both  sides  of  the  scattering 
surface.  Since  this  is  usually  an  indepedent  complicated  task,  differ¬ 
ent  approximate  methods  are  used.  This  work  has  employed  approximation 
of  the  physical  optics,  or  according  to  the  terminology  of  [2] ,  it  has 
assigned  a  "uniform"  component  of  the  current.  By  using  the  findings 
and  technique  [2]  ,  the  next  step  can  be  taken  to  compute  the  marginal 
waves  with,  regard  for  the  "nonuniform"  current  component. 

We  will  compute  the  diagram  of  the  scattering  field  from  the  ideally 
conducting  infinitely  thin  revolution  body  whose  surface  is  formed  by 
revolution  of  a  curve  of  the  second  order 

p=  P#:(l  +  tcostf)  (1)  _ 

during  its  revolution  around  the  axis  z  (fig.  1).  Here  e--eccentricity t 
Pq=(1+g ) f .  The  revolution  body  is  limited  by  the  plane  z=  -00’.  The 
phase  center  of  the  source  is  located  in  the  focal  point  of  the  curve 
of  the  second  order,  approaching  its  apex.  The  equations  which  link 
among  themselves  the  Cartesian  (x,y,z),  spherical  (R,0,$)  and  curvi¬ 
linear  coordinates  (u,  w,  v)  are  presented  in  [1].  The  mutual  arrange¬ 
ment  of  the  curves  in  the  section  of  the  studied  surface  with  plane 
y-0  is  shown  in  fig.  2a.  The  boundary  curves  are  characterized  by  the 
values  e=0  (circle),  e=l  (parabola)  and  c=®  (straight  line).  The  values 
c<l  define  the  family  of  ellipses;  e>l  defines  the  family  of  hyperbolas. 

As  the  analysis  made  below  shews  ,  the  obtained  solution  can  be 


Figure  1 . 

extended  to  the  case  where  the  source  is  located  in  the  distant  focus 
of  the  curve.  The  changes  that  are  needed  in  this  case  in  the  final 
expressions  are  indicated  in  the  appropriate  places.  Figure  2b  gives 
the  mutual  arrangement  of  th  curves  with  the  source  in  the  distant 
focal  point. 

We  will  obtain  an  explicit  appearance  of  the  equation  for  the 
revolution  body  with  curvilinear  coordinates  with  arbitrary  distance 
between  the  observation  point  and  the  focal  point  O'.  It  follows  from 
fig.  1:  pcos^«  -z  ,  in  relation  to  (1)  this  yields  p=Pq+ez. 

Publication  [1]  has  obtained  the  following  correlations  for 
arbitrary  z  and  p:  z=c[l+ch  u  cos  w)  cos0-sh  u  sin  u  sin  0]  and  p= 
c(chu+cos  w) .  By  using  the  link  between  z  and  p  by  curve  (3),  we 
obtain  the  following  appearance  of  the  equation  which  describes  the 
examined  surfaces  of  the  second  order: 


For  the  distant  zone  of  the  field,  the  studied  body  is  always 
located  only  in  the  second  or  third  quadrants  of  the  system  of  coordi¬ 
nates  z^O”h.  In  this  case,  in  the  initial  relationships,  the  parameter 
w  should  be  replaced  by  tt-w  for  the  second  quadrant  and  tt+w  for  the 
third  [1].  The  parameters  u  and  w  in  the  distant  zone  according  to 
the  modulus  are  much  smaller  than  1.  Accordingly,  we  replace  in  (2) 
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Figure  2a. 

l--straight  line;  2--parabola; 
3--circle;  4--ellipses;  5--hyper- 
bola 


Figure  2b. 

1- -straight  line;  2- -hyperbolas ;  3-- 
circle;  A--ellipses 

the  trigonometric  and  hyperbolic  functions  w-th  the  first  terms  of 
their  exponential  expansions  t  assuming  also  that  w  becomes  u+w.  The 
terms  of  a  higher  order  of  smallness  than  the  others  are  abandoned. 
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As  a  result  of  the  computations  we  obtain  an  equation  of  the  studied 
surfaces 


F(u,w,v)  s 2twtt>sin0$in  2uJ—  «‘(1+  tco$6)—  w»*(l — i cos 6)^.0,  (3> 

by  solving  it  in  relation  to  u  and  w,  we  find: 

«0  +  i  cos  0)  «  *  w  sin  0  sin  t*  + 

4.  |  2m| (1  +»cos0)—  —  »*-f  i*  sin* «  cos*  v)7 

u»(l  — ecosQ)  =  cMsinOsinv  + 
i  ]/2hJ  ( 1  —  *  cos  0)  —  u*  ( 1  —  s*  4-  sin1 0  cos zv)T  (5) 


The  square  root  sign  in  (5)  is  changed  in  the  transition  through  zero. 
The  root  is  made  equal  to  zero  at  the  point  of  contact  of  the  coordi¬ 
nate  ellipsoid  and  the  surface  of  the  body.  The  sign  before  the  root 
is  selected  from  the  condition  that  w  is  negative  in  the  second 
quadrant  of  the  coordinate  system  z^O' 'h  and  positive  in  the  third 
quadrant.  The  product  wsin  v  always  maintains  the  sign  "plus". 

The  differential  parameters  of  the  task  equal: 


»  =  — =[u(l  +  tcos0)— siwsln0sint>l:ltt>(l  —  scos0)  —  t  a  sin  0  sin  i»), 
du 

(6> 


q  «  __  (£UB>sin0cosu]:fa»(l  —  ecosO)  —  e « sin 0 sin  t>l.  (7> 


do 


From  (A)- (7)  we  find  the  necessary  functional  relationships: 


+  (1  + 


uw 


/-?_L  [ 2«| (1  — e cos 0)  —  «*(  1 — e* -f  **sln* 0 C03*vl 

u  ®  j "  i.  _ _ 


COSV 


V  »*— a|  I  a*  — 

«'(a*  +  ai) 


equal 


The  quantities  2"~va^ues  t^ie  coordinates  w  with 


t>  =  ± 


4-2- 


«,.;(!  --t-cos-0)  =  rwsin 0  ±]  2«o{l  —  scosH)  —  u* (1  t*). 

identical  to  the  value  w  on  the  edge  of  the  mirror:  ? 


The  quantity  B  is 


Ps  - 


2u 


2 

0 


(8> 


A 


1  4  «  cos  \|'0 


i.e..  (8)  is  the  equation  for  the  edge  in  curvilinear  coordinates. 

The  condition  2  yields  the  value  u  on  the  edge  of  the  mirror,  in 

the  section  of  it  with  plane  y-0.  By  solving  the  corresponding 
quadratic  equation ,  we  find 


cos* 


« 


ul  -  2 Ul 


0  l  +  t  COS 


(9) 


By  using  (8) ,  we  also  obtain 


u*  «=  2m* 


. . "  +  *r« 

sin*— -- 
0  i  +  *cos\f» 


The  "plus"  index  corresponds  to  the  upper  point  of  the  mirror  edge, 
the  "minus"  to  the  lower. 

The  integral  presentations  of  the  diagram  of  the  scattering  field 
[1]  for  the  given  case  has  a  different  appearance  depending  on  the 
combinations  of  integration  variables. 

The  variables  u,v: 


where 


at 


e-l*£'‘*  Vt(u)udu, 


(10) 


_ wdv _ 

V  (1  —  icos*i)  —  u* ( 1  —  i*-t  i*sin*iicos*ti) 


Variables  u,w: 


+ 

ikP±  i  G-ikeu'<p0(u)udu, 

nu0  u 


01) 


where 


<*>•(«)  = 


1 


.f 


wdw 


I  —  1  cos 6  J  a?)  o?— «'* 


Here  ~  OVs. »)  2pe<  p  ^ 


u*-H^"r*  M  ,u  --uniform  current  component 
on  the  surface  of  the  mirror,  i.e.  ,  TQ=2[nRQj  ,  n--unit  vector  of  the 
perpendicular  to  the  surface  of  the  mirror,  Hq- -unperturbed  value  of 
the  vector  of  magnetic  field  intensity  of  the  emitter  on  the  surface 
of  the  mirror; 
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As  the  emitter  we  will  examine  a  system  made  of  elementary  elec- 
trie  and  magnetic  dipoles  with  equal  moments  ,  oriented  by  the  x-nxis 
(electric  dipole)  and  the  y-axis  (magnetic  dipole) ,  This  emitter  is 
equivalent  to  the  .  flare  element  of  the  horn  emitter.  The  center  of 
the  emitted  spherical  wave  is  coincident  with  the  focal  point  O'.  In 
this  case 

£,«*•■  2 T  (9,  <p)  ( *'i  cos  T  4-  sin  <p) , 

7(0,  <r)  -  7„(0,  9)  (1  -cosO).  £0~  e°  1^- ;  (12) 

77, -  2T <tp.  g)  ?y  -  IM1 4-  cos i) sin Scos l  — 

— 7(cos*g  4- sin*  geos  sf)—fc  sinking ) ;  (13) 

Tq(6  ,<|>) --"scalar"  diagram  of  emitter.  The  vector  components 
perpendicular  to  the  curve  of  the  second  order  equal: 

_  ,  „  ,  .  I  +  C0S>1> 

a***  —  cosgcosx,  n„  ■» — sin^cos*.  rt,=»sin*.  tfix  ==»  -  — . 


x- -angle  compiled  by  the  tangent  to  the  curve  with  axis  z. 
After  simple  but  cumbersome  computations  we  find: 


/,  a  27  h|>,  g)  |f(sin  *  cos*  g  +  sin*  g  sin  (*  4  xfr))  — 

—  fsin2gsiny;Cos^x  4-  yj  +  Jcosgcos*  j  JL_I .  (14) 

6,  «=  27  (\|>,  g)  J|04-  s)cos<f  cos*  — (1  — *)  sin*  y  cos  (2| —  <p)  j  X 

XcosO  —  cos g sin  Osin  *).  05) 

»  —  27(t|>,  g)|(l  + '.) sin 9 cos* y- 4- (1  —  t)sin*y  sin(2g— <p)j. (16) 

The  formulas  linking  the  curvilinear  coordinates  with  the  angles  tj>  and 
g  are  produced  by  (8)  with  replacement  of  B  by  w  and  by  the  expression 


tgS- 


((W  _  u»)  cos  (i  4-  2ua>  sin  Q  sin  o)  sin  9  4-  cos  v  cos  9 
((»•  —  u»)  cos  U  4-  2«wsin  H  sin  c»)  cos  <p  —  2 uu<  cos  v sin  9 


07) 
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Analysis  shows  that  depending  on  the  position  of  the  observation 
point  in  space,  the  structure  of  the  scattered  field  changes.  For 


convenience  of  the  computations  it  is  expedient  to  divide  the  space  into 
a  number  of  regions  linked  to  the  geometry  of  the  mirror  and  the  posi¬ 
tion  of  the  observation  point.  Figures  3  and  4  show  the  division  of 
the  space  for  curves  of  hyperbolic  and  elliptical  types  respectively. 

In  each  of  the  isolated  spaces,  the  section  of  the  body  surface  with 
Fresnel’s  zones  has  a  different  nature.  The  analytical  presentation 
of  the  field  also  differs. 

Field  in  Illuminated  Region  of  Space 

In  this  region,  a  coordinate  ellipsoid  is  always  found  which 
passes  through  the  observation  point  and  touches  the  surface  of  the 
body  at  a  certain  point.  The  contact  condition  P~l  I  «  “  0> 


Att  »paeunt 
SoJIrttx 


Odna  it  pee  fan 


N 


ttUoJi™  \ 

..  —  . —  ■ —  _ \ 


I 


reoHempc  -  (’'•■’ni  necua  f> 
a  die  Kpaih.e  hump/ 


®  Kpaefgp  - r - i\ _ —  TTJZpUdfeS- 

n  Seana  #poM*u_  J[)  *“  y, - - 


Figure  3. 

I--Marginal  wave  of  edge;  II--Geometrical-optical 
and  two  marginal  waves;  III--Two  marginal  waves; 

IV--One  marginal  wave;  V--One  marginal  wave?  VI — 

Two  marginal  waves;  VII- -Marginal  wave  of  edge. 

In  addition,  at  the  contact  point  the  square  root  in  the  definition 
of  w  becomes  zero.  From  here  we  find 

cli</*  =  -j[±  i  '(s*  —  l)s  +  2«g !  (s4  —  1) cosO  +5a«j—  u\  J .  (18> 

,g%  „  “y1-  . 

1  —  t  ch  uk  cos  (i  jf 

The  "plus"  sign  corresponds  to  the  curves  of  the  hyperbolic  type,  the 
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’'minus"  to  the  elliptical.  The  value  c=l  is  boundary.  By  expanding 
the  square  root  in  (18)  into  a  series  and  by  limiting  ourselves  to  the 
first  terms,  we  find  the  following  ratio  which  is  correct  for  all  e: 


t  COS  fl  —  I  -f  — j — r 

«i  =  2u> - 1l+L 

*  0  **_!  +  .«* 


I  ‘“o'! 

^ ^  l,cas(i~1-fr+il0>-«+’"g) 

*  °(e cosfJ —  1)*  [  «*  — l-f  t  ^(H-cosO)]2 


Figure  4. 

I--Marginal  wave  of  edge;  II--Geometrical- 
optical  and  two  marginal  waves?  III--Two  margi¬ 
nal  waves;  V  --One  marginal  wave;  VI--Two  marginal 
waves;  VII--Marginal  wave  of  edge. 


The  scattering  field  with  values  e«l  (surface  of  the  body  is  close 

to  a  paraboloid)  in  the  narrow  region  of  angles  e  ~  0  must  be  analyzed 

separately.  This  is  associated  with  the  fact  that  with  small  values 

of  the  e  angle ,  the  entire  surface  of  the  mirror  lies  in  limits  of 

one-two  zones  of  Fresnel ,  therefore  the  contribution  to  the  field 

from  each  point  of  the  surface  of  the  mirror  is  significant.  Analysis 

2  2 

with  e  =1  was  made  in  publication  [8],  Accordingly,  the  case  f.  -1~Uq 
with  6~0  is  excluded  from  the  examination.  For  all  the  other  values 
of  the  angle  of  the  observation  point ,  the  obtained  results  are  correct 
with  e  =1  as  well.  Then  from  (8)  and  (19)  for  the  contact  point,  we 
obtain  the  following  relationships: 
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Figure  6. 


iff  A  -«-1tc  •  sin*  =  (t«-l)sinO_.  . 

g  2  *+!  5  ?  ’  5  **  «»+l~2eco8&' 

,  (t*  -4- 1'  ccs  0  -  2.- 

tos  **  -=  — - - . 

,(.  +  J  _  2t  cos  0 


The  angle  e'  (fig.  5,6)  is  linked  to  the  angle  of  observation  0 
by  the  correlation  sin  (O'  — 0)  --  u>  .  In  the  distant  zone  of 

the  field,  the  difference  between  these  angles  is  negligible. 


It  is  easy  to  obtain  from  the  geometry  of  the  task,  a  link  between 
the  angles  e'  and  ij< .  This  relationship  looks  like 


(20) 


i.e.,  is  identical  to  the  expression  obtained  above  for  the  contact 
point.  This  identity  characterizes  the  known  property  of  the  curves 
of  the  second  order  which  can  be  formulated  as  follows:  if  a  wave 
with  spherical  wave  front  is  emitted  from  one  focal  point  on  the  surface 
formed  by  a  curve  of  the  second  order,  then  after  reflection  from  this 
surface ,  a  new  wave  is  formed  also  with  spherical  front ,  but  with  center 
in  the  second  focal  point  (real  for  e<l  or  imaginary  for  e>l). 

In  order  to  write  the  integral  presentation  of  the  field  diagram 
we  will  examine  figs.  5  and  6.  The  right  side  of  fig.  5  and  6  shows 
the  pattern  of  curvilinear  coordinates  in  the  xz  plane  for  curves  of 
hyperbolic  and  elliptical  types  respectively.  On  the  left  side.it  is 
the  same  in  the  xy  plane.  As  is  apparent  from  the  figures,  for  values 
of  the  parameter  u  which  lie  in  the  intervals  uk~u+  (hyperbolic  curves) 
or  u+-uk  (elliptical  curves)  ,  the  intersection  line  of  the  arbitrary 
coordinate  ellipsoid  with  the  surface  of  the  mirror  is  a  closed  curve. 
Consequently,  the  Fresnel ' zones  in  this  case  are  ring-shaped.  The 
ellipse  with  the  same  value  of  u  intersects  the  x-axis  twice:  above 
the  point  of  contact  and  below  it.  According  to  this,  the  interval 
of  change  in  the  parameter  w  is  where  a=w  on  the  x-axis  with 

w>  wk,  and  ct2=w  on  the  x-axis  with  w<w^. 

On  the  remaining  part  of  the  mirror  ,  the  lines  u=const  are  not 
closed.  The  corresponding  Fresnel  zones  are  presented  in  the  form  of 
bands  which  rest  with  their  edges  on  the  edge  of  the  mirror.  The 
interval  of  change  in  the  parameter  w  with  fixed  u:  c^- B. 

The  integral  presentation  of  the  field  diagram  is  used  in  the 
form  of  (11).  In  the  first  interval,  the  limits  of  change  for  u: 
uk-u+  ,  for  W:  a1-a2.  In  this  case  0— ecosO) ^eusinOiK  (1—  e*)  («*  —  «*). 

In  the  second  interval,  u  changes  in  the  interval  u_-u+;  w- - in  the 
interval  c^-P.  As  a  consequence  of  the  symmetry  of  the  problem,  we 
will  do  the  integration  for  one  half  of  the  mirror  (for  example, 
the  right  in  fig.  5,6),  then  doubling  the  obtained  result.  We  will 
Introduce  the  normed  variables  of  integration.  In  the  integral  corres¬ 
ponding  to  the  Fresnel  zones  in  the  form  of  rings,  we  assume:  a) 
curves  of  the  hyperbolic  type  ( e > 1 ) : 
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i  ,  ,  ,  ...  2u*  (I  -4-  ‘  coi^0)  (t  cos 't  —  1)  —  n»j*| 

u  =  uk~  (  “1—  UD5  “  “ - - : -  ;  (21a) 

+  (1  4- >  cos  -»r*)  (-*  —  •+  *  «§) 

«■*  —  *|  -f  ( a*  —  a*)  p*. 

a.S*(.sln^-sin*±*)Y 


o  ft  u  sin  ^  ,  77 - 1 -  V  ?■  T'  « 

1- *> “  ,«nl 1  0 l: 


b)  curves  of  the  elliptical  type  (e<l): 


U*  ^u\-{  u\  -  “+)  **  1=1  (,  +  ,  COS  ^»)  (i  -  *  “o) 

xlo+.«.«o— 79»-“A 

,  /  -in  'ini?  +•  S  sin  — " —  )  ■ 


(21b) 


We  note  that  the  transition  from  curves  of  the;  hyperbolic  type 
to  curves  of  the  elliptical  type  in  the  obtained  relationships  is 
made  by  mutual  replacements  of  e-ljl-e  and  Therefore, 

further  analysis  is  only  made  for  hyperbolic  curves.  In  the  integral 
for  the  bands,  we  replace: 


:*  =  u\  +  ( ul-  u\)  t1  =  i  +  t  c--~  (a  +  bt*Y,  m*  =  a|)z?. 

(22) 

a  —  cos*--+2~0 ;  b  -  sin  G  sin  a  +  fr  =  cos*  . 


On  the  edge  of  the  mirror: 


B*  = - — - (aj-W*);  «»-  sin1——  ;  a, -s*-sin«  :  (23). 

y  1  + 1  cos  if,  2  2 

■  ;  COSMO- umtVtt  {.(0° 


■  ,8-** 
7=*r,in*  - IS  • 


A  sample  pattern  of  the  introduced  normed  coordinates  in  the  projection 
on  the  xy  plane  is  shown  in  fig.  7a, b.  Now  we  can  write: 

£?  «  L.  +  Li,  Lt  =  1  —  Yie"",*tM*  f  e1V,‘’  <t>(s)5rfs;  \24> 


where 


♦(*)« 


I 


I—  «  COJ  <t . 


%U)  -« 


I 


I  —  t  co*d 


■ 

I 


»•? 


_ £o^ _ 

l  1  ~Jr^* 


Y  -  2 kD  sin  0;  y,  —  y 
O  —  **»»!"  ^'o  . 

I  +  «  COS. 


(i»— l+«wj)b 


^--contributes  to  the  field  from  ring-shaped  Fresnel's  zones,  L^-* 
from  the  Fresnel  zone  in  the  form  of  bands. 

With  a  value  6-0^,  where  tgy  =  the  area  of  the  ring- 

shaped  Fresnel  zones  is  contracted  towards  zero,  and  becomes  zero. 
This  circumstance  determines  the  boundaries  of  action  of  relationships 
(24) :  0  j<  e  <_  0  .  The  interval  of  applicability  of  the  obtained 

formulas  will  further  be  presented  in  braces  ,  after  the  designation 
of  the  field  component,  for  example:  We  should  now 

compute  the  integrals  included  in  (24).  In  this  ease  it  is  expedient 
to  limit  ourselves  to  terms  on  the  order  as  compared  to  1,  since 

the  correctness  of  the  terms  of  higher  order  of  smallness  cannot  be 
guaranteed  by  the  initial  approximation  of  the  physical  optics. 

The  nature  6f  the  dependence  of  the  integrals  in  (24)  on  the  " 
large  parameters  of  the  task  is  determined  by  the  functions  Ms)  and 
*2  (tv,  >  therefore ,  one  should  first  clarify  the  nature  of  the  change 
in  these  functions  from  the  variables  s  and  t. 

We  will  clarify  the  nature  of  behavior  of  the  function  Ms)  de¬ 
pending  on  the  variable  s.  It  is  simple  to  show  th3t 


where 


COS0= 


V  ufi  —  a*  V  a*— 

»(aj— a^)(l  — tco*0) 


\  0  — £*T(  -  2s  Vi 


l—l **V  c*/ 


/]=•  V ** —  1  U*  sin*  S  (s  ccs  G  —  1)  (1  -f  e  cos  yk)  —  ( 1  —  e*4-  s*  sin*  S)a*s* — 

_ _  j_ 

— 2y<v(l—  2j**)e sin 0  (£i  —  1) K cos 0  —  l)  (1  +  tcos^*)  +  OjS*l|  *. 


It  follows  from  here  that  Ms)  is  the  function  of  the  combination  of 
variables  s/aT.  This  function  can  be  expanded  into  au'  exponential 


The  function  bQ  ^  was  written  above  in  an  explicit  form  [formulas 
(15)  and  (16)]  depending  on  the  variables  $  and  £.  These  variables 
are  associated  with  s  by  the  correlations: 

cosi-i-  -  _(t,~t)(K*  +  *c>)  +  2ua  .  cost  -  /Icoay-fisln?  , 

2  2t  («» +  u-t)  ’  s  )'A*  +  B*  ' 

A  —  “*  +  »*  — 


cos  (2j  —  qp)  = 
sin  (25  —  <p)  = 


B  —  utocoso; 

(>|«  —  fit)  cos  «r-  ~  sin  9 

( A *  —  fl«)  sin  9  4-2ABc«9_ 


At  the  contact  point  (s=0)  By  substituting  these  expressions 

into  $(0)  and  by  integrating  for  u  .we  find: 


Further , 


d  .  i  =  dfy.y  /a^  9u  ,  at  to  \ 
a*  *•*  **”°  dip*  \  du  ds  V  5»  <*»  \rax  ] 

dbhif  /dl  <*«  +  dcoup  \ 

*5  \  du  ds  V'Sj  dcosvds  (  n,  o 


+ 


But 


dw 


—  ~  0.  The  terms  which  are  proportional  to  ,  -  — 

d*V^  ,-o  *Ka»  Uo 


+t^90- (ecosQ-  .DO-V). 


are  also  equal  to  zero.  This 


is  associated  with  the  nullification  in  this  case  of  the  integral  for 
the  variable  y . 

By  fulfilling  the  operations  of  differentiation,  and  by  assuming 
in  the  final  expressions  s=0,  and  by  integrating  for  y,  we  find  as  a 
result 


2  -|/(i»-  1)(‘  cost)-  1) 

d)'(0)^  (7.4TiT^ru-2E  r  l~-f  *  cost.  X 

n[2^1SlS!(I  +  .cos^/,  #-27-(t*.  <r)0-«)0  +C0S^K)/f>l]. 
I  *<P  •  J 


The  intermediate  computations  which  are  necessary  to  obtain  the  des¬ 
cribed  relationship  are  omitted  here  because  of  their  elementary  nature. 

The  derivatives  of  higher  orders  can  be  computed  in  a  similar 
manner.  Further  we  should  substitute  (25)  into  the  first  integral  (24). 
The  integrals  for  s  are  expressed  through  known  functions  ,  and  each 
subsequent  result  of  integrating  series  (25)  is  an  order  smaller  than 
the  previous.  In  that  case  where  the  "large"  parameter  of  this  integral 
Yj  proportional  to  a2  approaches  zero  (which  occurs  on  the  boundary  of 
the  region  with  e-*e  )  ,  all  the  terms  in  series  (2i)  approach  zero, 
except  for  the  first.  These  features  of  the  obtained  expansion  $>(s)  show 
that  the  greatest  contribution  to  the  field  is  made  by  the  circle  of 
the  contact  point ,  and  consequently ,  the  expansion  into  a  series  near 
this  point  is  physically  justified.  In  this  case  there  is  a  possibility 
of  somewhat  simplifying  further  calculations  ,  taking  into  account  that 
the  boundary  line  u=u+  between  both  integrals  in  (24)  is  common,  and  the 
values  of  the  results  of  integration  coincide  on  it.  We  will  therefore 
take  each  of  the  integrals  in  (24)  once  by  parts.  The  values  on  the 
line  u=u+are  mutually  nullified.  We  obtain: 


£!.,  “  • 

®  (0)  e"'  -  -L  r'  *“*  f  e~lv'*'  * ; 

n  n  J  ds 

_ f  e-iv"*®i!»dr. 

*  i  * 


(26) 

(27) 


Here  E„°,  expresses  the  field  in  the  "geometro-optical"  sense,  since 
with  x-*0  in  (26)  the  first  term  does  not  disappear  and  is  precisely 
equal  to  the  "geometro-optical"f ield  [4].  As  will  be  shown  below, 

E*!we  corresponds  in  nature  of  the  field  to  the  ''marginal"  waves. 

r  *V 

Based  on  the  obtained  correlations  ,  it  is  easy  to  compute  the  value 
Ej°  .  After  integration  for  s  with  regard  for  the  equality 

|  e"“*  ds-|/^lC(Yi)  +  i5(Y,)l. 

where  C(x)  and  S(x)--corresponding  Fresnel  integrals  [3],  we  obtain 
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{0  *r  0f)  «  (£,  -j-  Et)  e  *  /,  f  +  C-'H!  ik  '  , 


(28) 


where 


£. 


dip 


«ino 

h3*  tC(Yi)  +  1  S  (V,)| 

Vn  W>.!(tI+l)«ao~2  )» 


In  addition , 


F--  ,  r  rn.  <‘-D»  W-l) (»cosu- 1)  C(Yi>+‘S<Y«> 

£«.t ~~7* f )  0 - ((tt  f  iTTmo-^T* - v  -^  •  ( 


Formula  (28)  shows  that  the  geometro-optical  v  E^  whose  ampli¬ 
tude  does  not  have  a  frequency  relationship,  can  be  viewed  as  the 
first  approximation.  The  amplitude  of  this  wave  repeats  the  nature  of 
the  diagram  of  the  emitter  depending  on  the  angles  and  <)>.  The  compo¬ 
nent  E2  which  has  a  frequency  relationship  is  a  correction  for  the 
geometro-optical  field.  This  correction  is  also  linked  to  the  contact 
point.  Its  amplitude  depends  on  the  degree  of  nonsymmetry  of  the 
emitter  diagram  for  the  angle  <J>.  The  diagram  of  the  emitter  may  also 
be  nonaxisymmetrical  since  until  now  this  circumstance  was  not  a  re¬ 
striction  on  the  analysis.  The  correcting  wave  has  fche  contact  point 
as  its  "source." 

The  quantity  EQ  characterizes  the  cross-polarization  component 
of  the  field  which  is  also  associated  wit!  reflection  of  rays  at 
the  contact  point.  This  component  has  a  frequency  dependence  and 
becomes  zero  at  8=0  (as  a  consequence  of  symmetry  of  the  mirror)  and 
with  9=er  (in  light  of  the  nullification  of  L^) .  With  e=l  (paraboloid) 
the  cross-correlation  component  is  also  missing  (the  current  line  does 
not  have  a  y-component) . 

To  completely  characterize  the  field  in  the  region  0-8r  we  should 
also  examine  the  component  (27)  since  precisely  this  component  must 
"extinguish"  the  action  of  the  geometro-optical  wave  on  the  boundary 
of  the  region. 

In  a  similar  manner  we  first  determine  the  explicit  form  of  4s  (t). 
We  expand  the  function  b.  .  into  an  exponential  Taylor  Series  for  z 
near  z=l  (fig.  7,  edge  of  mirror)  ,  after  which  we  integrate  the 
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obtained  series  for  z .  We  find 

-  r-T^r,[ w  arc  sin  ~ £ v,U  +  •  •  ■]• 

Here  in  order  to  simplify  the  analysis,  we  introduced  the  restriction 
that  bQ‘  does  not  depend  on  z.  This  is  equivalent  to  the  hypothesis 
that  the  diagram  of  the  emitt  r  T(\|»ty)  is  axisymmetrical ,  i.e.  ,  depends 
only  on  the  angle  <ji .  According  to  (17)  we  have  on  the  edge  of  the 
mirror: 


tgs~ 


( 1  —  2**)  sin  y  4-  2/  \‘  T—  <*  cos  if 
( 1  —  20)  cos  ip  —  2 1 )  1—  i*  sin  if 


e-KlS-f)  _  ( j  _  gji  +  8(«  4.  { 4*  J  ’  I—  }*)  e,f • 

By  substituting  this  definition  of  $  and  into  (15)  and  (16)  t  we 

obtain 

6jp  (I)  =  IT  (if0,  ()  jeos (f  [a  cos  6  4-  cos  (0  +  \f#)  -f  8(1  ~  4)  sin*  X 
>:  cos  8  /*  ( 1  — /*)  +  2/*  sinOsin  \j:#j  4-  sin  <f  |sin  6  sin  2/J  1—  (*  4- 
4-  (I  — «)  sin* -y  cos 04/ ]  T^7*(l  —  2/*)]} ; 

(/)  =  —  2T(%,f)  |sin  4-s cos  to— 0  —  s»n*  y  8<s0  —  **)]  + 

.e  4  (1  —Osin*  %  t  |/|  — l*(l  —  2/*)cos<r  j .  (#>) 

4.  J 


Since 


1  *0  ( °?  “  cf)(i  4-  ecn5  4»)  1  1  + :  cos to 


Zu’ib  L  sin*« 

2(c,  —  bt*)  r  :  cos'll  Mb,' 


l/T+W  2(*,-bfi)  4- =  coo  to  *K%  1 

—  1  y-^TTT  “cos*  4,  4-  V  4-  :>rp  (o,  —  bn)  cosy  (77  Oi  Jkp' 


4  i- 

m sl 


then  one  can  find  the  following  term  of  the  expansion  4>j,(t)  etc.  in  the 
explicit  form.  The  quantity  Xq  that  is  included  here  equals: 

.  |'o4  W«1-g«  -i-»*i«4-  1  Osip.  ti>(:  -  C»S.I)  . 

*  ~  " ~  ~  ''“ZZSZ  ~  "  * 

2  1  o+  bt’  |  o, +  {>';« 

*'=»<**—  1)*;  x,(i)«o. 
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Now  (26)  adopts  the  appearance 


—  I  .♦  <1 

05  - -Tibsij  [*:r.w»«v 

jrf,. 


It  is  sufficient  to  be  limited  only  to  the  first  term  in  the  brackets. 
The  following  terms  of  the  expansion  provide  a  correction  for  the 
negigible  order  of  smallness,  of  which  we  can  be  convinced  by  direct 
integration  by  parts.  For  similar  considerations  we  should  abandon 
the  terra  ^Larcsin.v.  (we  recall  that  3/3t  must  be  fulfilled  with 

dt  # 

regard  for  the  dependence  of  b.  on  <Kt)  and  £(t),  direct  application 
of  this  operator  to  b0F^  is  unlawful).  By  differentiating,  we  also 
find 


.  d 

—  arc  sin  xa  — 
it 


sin  to  ( 1  —  s  cos  it)  /  £  yf 
y  j__|*  \a-\-bl* 


•1*  1  o.  j 

n*  t  b’i:  )  ‘ 


(30) 


In  this  case, 

A  —  ]  a  4- b  :  [l'o  ( a .  +  b ')  —  )  o.,(o --  />))  ; 
A"  —  Vot+  b' :  (}  a  (a%  ’4-  b')  —  l  ws(a-f  b)]  - 


With  regard  for  what  has  been  said ,  we  obtain 


£»?***  —  sin  ^0e 

T  n 


—I  *£«' 


r+  (■  -m* 

)  c  \a+  bt 


A "  l  a,' 


dt 


6/s  uj-fc'/s'i  i—  /* 


With  ^ i  the  terms  of  the  required  order  are  obtained  when  only  the 
first  term  is  kept  in  b^P^ 

where 


=  £  cosO  +  cos (6  +  \J?0);  A®  =  1  +  ? cos \f0. 


In  the  portion  of  the  region  of  the  observation  point  angles,  in 
which  the  parameter  y>>1  we  obtain 


2  —I  te«<* 


__  V  "Vw  t  ^ 
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-J 


! 

\\ 

i. 

i  ' 

r 

t 

1  i 


’*1 

y  ; 

V 

fci 


|_  tt  m^TTIy.1.  •  _?  -'-- 


Here 


\ 

(a,b)  **  j*< 


-i >fi« 


it 


(«  +  M*)r  l~f» 


We  use  the  obvious  correlations  »~m> 

«+bt* 


?*T  XJ 


•  -I  ij?  )•**»> 
e  *  <fr. 


By  substituting  it  into  the  original  integral  and  by  integrating  for 
t  [3]  we  find 


I(a,bh 


aj> 


V 


This  integral  is  expressed  through  the  cylindrical  functions  from  two 
actual  variables  (Lommel  function)  Uq(c,x)  and  (c  ,x)--f irst  and 
second  order  respectively  [5].  These  functions  have  been  tabulated 
and  well  studied  (6,73.  We  can  write: 


/(«,*)* 


iu 


|K  a(°+ 


|r[£/.(c.f .  f )  +  w,(«. -f  •  f )-TJ*(i)l: 

(Va  +  fe-Vrfl>1 


/!,  *s 


(32) 


By  using  the  asymptotic  presentation  of  Lommel's  functions  in  the  form 
given  in  [7] ,  we  obtain  in  our  designations 


x[~— c(^)  j-  +  is(5 


b 


For  7  >  >1 


/(a.  b)t 


2»«Y 


^i+i£L.-+]. 

c  +  *  x  ya  J 


i’ 
>•  ' 


h  i 

t 


where 


f(Y)“rhll|D"l+i'T,'26'YIc,(“)_ 

—  iSt^)]siSn^  *»• 

Ci(x)  «*  ~ — c (*):  -j-"5 {x)' 


Now  for  the  unknown  component  of  the  field,  we  can  write: 
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where 


-I'K+i 


+  /\+>  e 


-l  *»<.+-•  T  | 


I-  a  4*  6  )  fl|  +  6 

The  phase  multipliers  in  the  first  and  second  terms  (33)  show  that 
the  ’’source"  of  the  first  wave  should  be  considered  the  edge  of  the 
mirror  -D,  since  the  coordinate  ellipsoid  which  is  determined  by  the 
parameter  u_  passes  through  it ,  while  the  "source"  for  the  second  wave 
is  the  edge  >D  for  a  similar  reason.  Thus,  (33)  characterizes  the 
known  marginal  waves  in  approximation  of  the  physical  optics.  The 
amplitudes  of  these  waves ,  in  contrast  to  other  asymptotic  methods , 
are  continuous  everywhere  and  do  not  have  peculiarities  near  caustic  curve 
surfaces.  The  nature  of  the  change  in  the  diagram  of  the  field  near 
the  special  points  is  determined  by  the  F  functions.  By  using  the 
asymptotics  of  Fresnel's  integrals  with  large  values  of  the  argument , 
we  obtain 

'Pr)=7T  ** 

'(-Pb-iir  -"ith 

Accordingly ,  in  that  part  of  the  examined  region  of  angles  of  the 
observation  point  where  the  indicated  parameters  are  fairly  great , 
we  have 

p»>  -  rk  •  (34> 


Below  (see  fig.  15)  is  the  course  for  the  actual  and  imaginary  parts 
of  the  F  function  and  its  asymptotic  presentation  for  a  certain  critical 
case.  It  is  apparent  from  the  figure  that  the  narrow  region  of  angles 
6  which  adjoin  e  is  transitional.  In  this  transitional  region  of 
angles ,  the  marginal  wave  changes  the  order  of  dependence'  on  the  large 
parameter.  The  ’’flash"  of  the  marginal  wave  that  is  observed  in  this 
case  partially  "extinguishes"  the  geometrical-optical  wave  near  the 
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boundary  of  the  region. 

By  uniting  the  obtained  results,  we  can  write  for  the  scattered 
field  with  0<  o 


s- 


,,  <♦»  -  l)v>  (<  COS  1)S/*  (Ct  (VO ±  1 5t  (v,)l  t  c“l  *4  /t  + 


+ 


*.U~> 


I(t»+  1)coj8  4-  2eJ* 

l 

—  I  *s« 

e 


(35) 


The  cross-polarization  component  of  the  field  provides  expression 
(28a).  The  complete  field  is  obtained  by  adding  to  (35)  the  field  of 
the  emitter.  The  scattered  field  is  the  sum  of  the  geometrical-optical 
wave  and  the  two  marginal  waves,  Althought  the  marginal  waves  which 
are  computed  only  for  the  uniform  current  component  generally  inac¬ 
curately  describe  the  actual  picture  of  scattering,  neverthless ,  in 
the  examined  region  of  observation  point  angles  in  which  a  geometri¬ 
cal-optical  *  field  exists  which  is  dominant  in  value,  (35)  apparently 
yields  a  satisfactory  approximation  (naturally,  for  kD>>l).  In  the 
direction  of  the  boundary  of  the  region  (e=0r),  it  follows  from  (33) 

i  — i*«J  , 

e-rpmi  ifj  )  s  —  T  (ih  ) - - — - L  e  + 

t'r I  '  ,*4- l_  2tcosOr 

sin  I"  1  -M  cos  e“'  *ta- +  1 T  4. 
j'ny  sin  Or  [  s  +cos^,  l+icos'r. 

It  is  obvious  from  here  that  the  marginal  wave  extinguishes  the  geo¬ 
metro-optical. 

The  obtained  relationship  (35)  is  incorrect  near  0~O  since  when 
the  expression  is  derived  for  the  marginal  waves,  it  was  assumed  ~  >1. 
With  0  close  to  zero,  the  marginal  waves  are  transformed  into  toroidal 
waves  of  the  edge  of  the  mirror  which  in  this  case  is  excited  by  the 
source  in-phase.  To  clarify  the  characteristics  of  the  wave  of  the 
edge  in  the  integral  presentation  of  the  field  diagram  (31),  all  the 
terms  in  bQ  have  to  be  taken  into  account.  The  integrals  obtained 
in  this  case  are  known  [3].  By  omitting  the  intermediate  computations 
and  by  only  entering  those  terms  which  are  significant  even  with  e-l, 


*  " 

-ttnfl  -I  - 

e  +  * 


' 

m 


(36) 
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we  can  write  the  following  expressions  for  the  wave  diagram  which  is 
"emitted"  by  the  mirror  edge: 

i'C" P  < 4)  — ; [''«•'•(!) '  *5- 1  “* (f)- 

(37) 


sin  **’ 

£;■”  (»« «)  -  n«  — - 


V  Y  Y  ?  v 

s,n  '2'  2  cos  2  -<  *«■+“'**  ,oos 

- - - j - e  ,  (J8) 


(i/ 


where 


£0  =  6JT[)  aas— (s -f  cosi|)0)  sinOj ;  =  (1  —  :)sin--y  cos&; 

*7  =  — (1  — 2)sin2^- ;  b%  =  sin  0  p.  +  cost):0  -f  fl«  +  (t*Z  l)  ■  -  ]  . 


Here  Aq(x)  and  A-^(x)  --lambda-functions  of  Bessel  [3].  In  computing 
the  field  near  e~0  ,  the  marginal  waves  in  (35)  must  be  replaced  by  a 
wave  which  is  defined  by  expression  (37).  A  similar  component  (38) 
must  be  added  to  the  cross-correlation  component  (28a).  It  charac¬ 
terizes  the  action  of  the  mirror  edge.  We  note  that  with  e+1  ,  formula 
(37)  reflects  the  effect  of  amplification  which  is  inherent  to  the 
parabolic  mirror. 

In  this  case  we  obtain 

2tg—  _l*c  2  _  |  1 

S"T.  ,39) 

As  already  indicated,  the  case  e  =1 ,  in  light  of  its  practical  importance, 
was  analyzed  separately. 

We  will  .  now  write  formulas  which  determine  the  field  diagram  in  the 
remaining  space.  In  light  of  the  fact  that  the  nature  of  the  compu¬ 
tations  is  similar  to  that  already  given,  and  in  this  area  of  space, 
the  field  is  only  created  by  marginal  waves  which  in  the  approximation 
of  physical  optics  can  be  viewed  only  as  auxiliary  material  for  appro¬ 
ximate  calculations,  the  calculation  details  are  given  in  the  appendix. 

For  the  part  of  the  illuminated  region  in  space  that  is  included 
between  the  boundary  of  the  geometrical-optical  field  (6=0  )  and 
the  line  of  the  boundary  of  light  and  darkness  ( G  — n - ^ q )  .  we  find  the 


following  analytical  presentation  of  the  diagram  of  the  scattered 
field  [formula  (A. 5)]: 


1'o>=~Oo[/ 


-<  *f«L  + 1  -- 

a  * 


where  P^_^  and  P^  are  defined  as  already  known  functions. 

In  the  middle  part  of  the  region  where  the  arguments  of  both  F 
functions  are  fairly  great,  according  to  (34)  we  have  P(4.,«s  y^' 

In  this  interval  of  angles,  shading  of  one  edge  of  the  mirror  by 
the  other  is  possible.  The  changes  that  develop  in  this  case  in  the 
diagram  of  the  field  are  not  taken  into  account  in  the  approximation 
of  the  physical  optics.  However,  one  can  draw  definite  conclusions 
from  the  results  of  a  series  of  numerical  calculations  which  are 


presented  below. 

It  is  easy  to  verify  that  substitution  of  the  values  0=0r  in  (40) 
turns  this  formula  into  (36). 

The  obtained  relationships  (35)  ,  (37)  and  (40)  completely  charac¬ 
terize  the  diagram  of  the  scattered  field  in  the  adopted  approximation 


for  the  part  of  space  illuminated  by  the  sources, 
boundary  of  light  and  darkness  we  obtain:  (r—j0), 


On  the  line  of  the 


Field  in  the  Shaded  Region  of  Space  (”— 'to^Os^n). 

In  the  shaded  region,  the  characteristic  division  of  the  mirror 
surface  into  Fresnel's  zones  has  the  appearance  shown  in  fig.  8. 

The  line  u=0  intersects  the  surface  of  the  scattering  body.  Part  of 
Fresnel's  zones  is  ring-shaped,  part  has  the  appearance  of  a  band.  The 
picture  of  Fresnel's  zones  on  the  surface  of  the  body  is  similar  to 
the  examined  case  0_<  0  <_  0r.  We  therefore  make  an  analysis  in  the 
same  way.  The  integral  presentation  of  the  scattered  field  is  written 
in  the  form  a 

£?  =  _  i  — 2-  \  C~ikCU'  («)  udu  +  J 

*  nuo  j 

.  +  “ 


Here  in  H'q(u)  the  limits  of  integration  for  v  are  constant:  “7  “  7  ♦ 
We  will  take  the  entered  integrals  once  by  parts.  As  shown  in  [1], 
the  field  with  u=0  equals  the  field  that  is  created  by  the  emitter  at 
the  observation  point,  but  with  the  opposite  sign.  The  values  of  the 
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Figure  8. 

results  for  integration  by  parts  for  both  integrals  on  the  line  u=u+ 
are  mutually  nullified.  The  integral  presentation  of  the  diagram  is 
written  in  the  variables  t  and  z  which  are  introduced  according  to 
formulas  (22)  with  mutual  replacements  of  and  02: 


JV-* 


-lto4  1 


(41) 


u> » **+(P*--3tI,s*;  ®*(0  =  i_e»Tof  1  T^T  ' 


The  replacement  of  and  a  2  is  caused  by  the  fact  that  the  upper  and 
lower  edge  of  the  mirror  now  fall  in  the  third  and  second  quadrants 
of  the  coordinate  system  ho'^  respectively,  and  not  in  the  second 
and  third,  as  occurred  in  the  illuminated  region.  The  parameter 


jc*  = 


Va  -f  bt*  Vat  +  b't*  —  )'oa,  +  <«(:-  cos .))  sin  \}-<, 


(42) 


2  Y 0+  bt1  f  a,  +  b'/» 

By  analogy  with  the  previous  calculations,  we  can  write: 


BnoJta  =  FO  *  £0  • 


f 


f  e-iV,‘  bZ  4  *rc  sin  A'» dt  +  0  (t) 

(1—  C  cosu)*'  di  '  ^  ; 


i-ncwiM 
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Haiti)* 


d  .n  $ 

It  is  easy  to  obtain  sin  x,  =  cos  ,\q  ,  and  consequently  '■£  arcsm.v,  ---  arc  sin  .vn 

It  follows  from  here  that  the  final  results  can  be  written  in  the  form 
of  (33),  but  with  a  change  in  the  sign  in  front  of  the  formula: 


Hr  {(*-*«) 


'}  a0[p(-> 


-I  kcu_  +  I  —•  „  -i  kcu j. 

e  4  +P{«  e 


j.  (43) 


Formula  (A3)  is  inapplicable  directly  near  the  line  of  the  shading 
axis  0  =ti  since  the  hypothesis  y>>1  is  violated  which  is  the  basis  for 
its  derivation.  In  this  case  we  should  act  in  the  same  way  as  in  the 
derivation  of  formula  (37),  i.e.,  take  into  consideration  the  main 
terms  bQ  .  By  abandoning  those  terms  wliicl^  are  not  important  with 
0~tt,  we  find 


*sr  — r  <*>  T?tr  '<■' 


-!  *fu+  -  I  -j 


04) 


Expression  (44)  characterizes  the  toroidal  marginal  wave  of  the 
in-phase  excited  edge.  The  field  amplitude  of  this  wave  is  an  order 
higher  than  the  amplitude  of  the  normal  marginal  wave.  Consequently, 
in  a  strictly  "backwards"  direction,  the  field  diagram  has  a  maximum 
which  is  comparable ,  and  in  other  cases ,  larger  in  size ,  than  the 
field  amplitude  during  emission  "forwards."  This  is  the  so-called 

fn 

Poisson's  loop  in  the  field  diagram. 

Extension  of  the  Solution  Results  to  the  Case  of  the  Source  Located 
in  the  Distant  Focus 

The  geometry  of  the  problem  and  the  employed  coordinate  systems 
when  the  source  is  located  in  the  distant  focus  O'  of  the  curve  generatrix 
are  shown  in  fig.  9  (as  an  example,  the  generatrix  is  taken  in  the 
form  of  a  hyperbola) .  The  equation  for  the  surface  of  the  scattering 
body  in  this  case  is  described  in  the  form  4 


P  “  : ( 1  — J cos *4"' );  p;»(l  j^/.  (45) 

The  center  of  the  cylindrical  coordinate  system  is  taken  at  the 
point  O',  the  angle  if* '  is  counted  from  the  z-axis.  The  parameter  e 
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Figure  9. 

--eccentricity  of  the  curve;  f --distance  between  the  focus  and  the 
nearest  apex  of  the  curve;  F--dis ranee  between  the  focus  and  the  distant 
apex  of  the  curve.  We  recall  that  the  mutual  arrangement  of  the  curves 
in  this  case  is  shown  in  fig.  2b.  By  comparing  (1)  and  (45)  ,  we  note 
their  identity  with  replacements  of  f+F  ,  c  j  -e.  The  angle  <(> '  is  the 
angle  at  which  the  current  point  N  of  the  body  surface  is  visible  from 
the  distant  focus  of  the  curve ,  at  angle  e  ,  the  angle  is  elevated  at 
which  the  observation  point  M  is  visible  from  this  focal  point.  The 
angles  and  f  are  associated  by  the  relationship 

e  2  e  -r  I  *  2  ,  , 

Figure  10  presents  the  designations? of  the  employed  parameters  of 
elliptical  and  hyperbolic  curves  wi  th  arrangement  of  the  source  in  one 
of  the  foci.  Comparison  of  the  nature  of  division  of  the  body  surface 
into  Fresnel  zones  in  both  cases  ,  both  for  the  illuminated  part  of  the 
space  (fig.  5  and  correspondingly  cig.  11)  p  and  for  the  shaded  (fig.  8  and 
correspondingly  fig.  12)  ,  shows  their  basic  identity. 

Taking  into  consideration  that  the  listed  factors  are  definitive 
in  deriving  formulas  for  the  field  diagram,  one  can  draw  the  following 
conclusion:  all  the  relationships  obtained  in  the  previous  sections 
maintain  their  meaning  and  transmit  the  characteristics  of  the  diagram 
of  the  scattering  field  when  the  body  surface  is  irradiated  from  the 
distant  focus  of  the  curve  with  the  replacements: 
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Figure  10.  l--Hyperbola i  2--Ellipse; 
3- -Asymptotes . 

t~Fi  e«- - £;  Po~~P'0  (46) 


and  with  counting  of  the  angles  t|»  and  6  from  the  distant  focal  point 
of  the  Curve.  We  recall  that  the  previously  obtained  formulas  were 
written  for  curves  of  the  hyperbolic  type.  In  the  transition  to  curves 
of  the  elliptical  type  ,  as  before  ,  additional  permutations  should  be 
made  in  them' 

c  -J-  1  — ►  1 — s  at-*-at. 


For  economy  of  space,  the  corresponding  finite  expressions  are  not 
written. 

Certain  Solution  Results 

Analysis  of  the  solution  shows  that  in  different  regions  of  space, 
the  problem  has  different  parameter  characteristics.  The  quantities 
y  ,y /b  and  ya2  can  be  isolated  as  these  parameters  with  "large”  para- 

vT 

meter  of  the  problem  kPQ.  On  the  boundaries  between  the  isolated 
regions  (0;  0r  it  - q  j  -n )  one  of  these  parameters  is  nullified.  The  order 
of  the  dependence  of  the  field  amplitude  on  the  main  large  parameter 
of  the  problem  IcPq  changes  correspondingly.  The  boundary  lines  of  the 
regions  correspond  to  the  caustic  curves  of  the  field.  At  a  sufficient 
distance  from  the  boundary  lines  ,  all  the  characteristic  parameters 
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Figure  12. 


are  "large." 

The  findings  confirm  the  principle  of  locality:  the  scattering 

field  has  that  phase  structure  as  if  it  had  been  created  by  locally 
reflected  rays  of  the  source  from  the  corresponding  points  on  the 
mirror  surface.  The  directions  near  the  axis  of  body  symmetry  are  an 
exception. 
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Figure  13. 

1 - - db ;  2- -marginal;  3- -region 

,  f„,ned  svstea  of  formulas  makes  it  possible  to  compute  the 
The  obtained  system  including  the  transitional  regions, 

scatterrnc  field  rn  lraation  0f  physical  optics.  The 

» ... ••  >"*  ,oi  r.nS 

mate  calculations  of  the  scattering  field  from  hyperbolic  and  P 

tieal  auxil  ary  mirrors  used  in  standard  two-mirror  antennas ,  as  we 

thr  radiation  field  from  the  main  parabolic  mirr  . 
as  to  compute  the  r  diagram  were  made  for  a  number 

Calculations  of  g  integrals  were  computed 

of  examples.  In  publicatio  .  .  13  and  u  present 

by  the  methods  of  numerical  integration.  g 
y  *  fb*  matter  ins  field  which  is  computed  fiom  tne 

rri—Lf^rrhy^erbolic  reflector,  "« 

analogous  curves  IM  ->»•* •  J  the  planes  E(»- 

^HU-T’ln  the “  plane.  a  dip  appears  in  the  diagram  which 
could  be  e^lalntd  by  the  corresponding  orientation  of  current  e  emen^^ 

flowing  on  the  edge  of  «h.  - £>£ ,  d  ‘ thl ’^tion  point, 

in  which  the  current  vector  is  directed  towar  ^  kl 

the  current  flowing  perpendicularly  to  the  edge  era 
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l--ellipse;  2--hyperbola. 

lies  beyond  the  limitis  of  accuracy  which  can  be  guaranteed  by  the 

approximation  of  physical  optics.  Figures  13  and  1A  also  give  different 

components  of  the  field  which  were  computed  separately  from  formulas 

(28),  (33),  (37)  and  (AO),  This  separation  of  the  contribution  of 

each  component  makes  it  possible  to  orient  oneself  in  the  methods  of 

controlling  undesirable  radiation  in  the  selected  directions.  These 

same  figures  illustrate  by  the  dotted  line  (in  the  region  Ilia)  the 

field  from  one,  unshaded  edge  of  the  mirror.  As  is  apparent,  the 

difference  between  the  interference  field  of  the  "light"  edge  and  the 

shaded,  as  compared  to  the  field  of  one  "light"  edge,  is  not  great,  and 

may  not  be  taken  into  consideration  in  practical  orienting  calculations. 

One  can  draw  the  conclusion  from  here  that  although  in  the  approximation 

of  physical  optics  in  the  field, rays  are  present  which  "puncture"  the 

mirror  metal ,  their  amplitude  is  significantly  lower  than  tne  amplitude 

of  the  actually  "shining"  edge.  .  ;* 

Figure  15  presents  the  nature  cf  the  change  in  the  furction 

for  the  examined  example.  The  dotted  line  shows  the  value  of  the 

analog  of  this  function  of  the  quantity  A  which  follows  from  the 

!  I  a  l 

method  of  the  stationary  phase.  Kith  #  0r|— 'J°* 
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Figure  16  presents  examples  of  computing  the  field  diagram  for 
a  number  of  specific  values  of  parameters  c  and  kD.  It  is  apparent 
that  the  "conjugated”  ellipse  and  hyperbola  (••  »-r  —  Ip  have  practically 
the  same  scattering  diagrams. 


Conclusions 

A  system  of  formulas  has  been  obtained  which  makes  it  possible 
in  the  approximation  of  physical  optics  to  find  at  any  point  in  space 
the  characteristics  of  the  scattering  field  of  a  spherical  wave  on 
finite  revolution  bodies  formed  by  curves  of  the  second  order.  The 
found  solution  can  be  used  to  make  the  next  step  in  computing  the 
nonuniform  component  of  the  current. 

The  analysis  made  shows  the  effectiveness  of  the  method  of 
physical  coordinates  in  solving  scattering  problems. 

The  results  of  the  work  may  be  used  in  different  areas  of  antenna 
engineering,  for  example,  in  computing  the  scattering  field  from  mirror 
antennas  or  re- emitters  of  the  appropriate  shape. 


Appendix  1 . 

We  will  clarify  the  nature  of  the  relationships  that  determine  the 
field  in  the  region  of  angles  er<_  e  <_  tt-Vq.  In  this  region,  there  are 
no  rays  that  are  reflected  from  the  mirror  according  to  the  laws  of 
geometric  optics.  The  integral  presentation  of  the  field  is  written 
in  the  variables  u  and  v.  Expression  (lO)with  replacement  of  u  and  v 
by  normed  variables  t  and  t  (fig.  7b)  according  to  the  formulas 

a*  «*!$.  +  ( «!—«+) <*  u  ~  4-[v(/)  — j-]*. 

adopts  the  appearance 

9  -I  ktu*  ' 

— fe 
4 

where 


l  (0  =  arc  sin  (£0  (/)  ]  1  —  t*  1 1 ;  sin  £  (/)  —  cos  v  (f); 

k(t)~z  sinQ-i-i±^l-  ; 

I  \  +  b7* 

'  ..  1 
q(t)  ~  \  \,12U5(1— icosO)— u*(l— :*-j-i2sins0sin5^T))  *  u\{r. 
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We  assume  //  ^(/hin-(/)r,  then 
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(A.  2) 


Similarily  to  the  analysis  given  above  ,  we  expand  b.  into  a 

U  t9 


Taylor  series  for  the  exponents  y  near  the  edge  of  the  mirror  ini im] 

It  follows  from  the  previously  obtained  correlations 


— i* 


dbi, 


dy  '"**  *  (0 1  (1  -  2/*)  (/)  4  t  ( I  •  is)d-rr^  |  **  i*n 
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After  integration  for  y  we  determine 


.9(0  -  r; 


i 


|(0(l -tcost*)l 


j  [arc  sin  1* (()  sin  5  (/)  — 


[l  -V'l— **(0*in*S(0-*(0(l-cos5)J  + 

W*  Kp 


The  numerical  presentation  for  q(t)  should  be  substituted  in  the 
original  integral.  It  is  easy  to  be  convinced  that  the  second  term 
of  the  expansion  already  yields  quantities  of  a  negligible  order  of 
smallness.  By  integrating  the  remaining  expression  once  by  parts  we 
obtain 


-I  k(U*i 


■  j  e_l  yr  (0 17  larc  sin  k  (0  * (0  ■ 


(A. 3) 


n(J  — scosd) 

—  arc  sins  (01  <0  +  O^-j 


By  differentiation  we  find: 


a  ...  itt- MVa.±?+Vai-  b  YW)  . 

-  arc  sin  ft  (0  sin  5(0  -  “  ‘  *  —^Zbt%)  {a%  +  b't .)  |  'P7« 


dt 


a  .  .  i  an.  —  i»d  o  t>  v'fii  —  j  x  Ijgfij 

'  •  arc  sin  * to  -=  -  s,n  (a  +  61*)  V  CT 


In  addition , 
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By  substituting  these  correlations  into  (A. 3)  we  arrive  at  integrals 
which  coincide  in  form  with  those  previously  computed  [see  (32)  and 
so  forth].  One  can  therefore  immediately  write  the  result  of  inte¬ 
gration.  Since  y>>1  always  in  the  examined  case,  then,  by  using  the 
asymptotics  of  the  Lommel  function ,  we  find 
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Passing  in  this  expression  from  the  letter  designations  to  the 
corresponding  functions  of  the  angles  ,  after  simple  computations  we 
obtain 


(*— ?o)|  -=  — Oo 


—I  koj^_-H 
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4  +  P(+) 


-|  KCM^-I  " 

e  +  < 


]• 


(A. 5) 


which  is  identical  to  (33)  with  accuracy  to  the  sign. 

Consequently,  the  marginal  waves  in  the  region  of  angles  0 < © < it 
are  described  by  the  same  analytical  relationship.  In  this  case 
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BROAD-DIRECTIONAL  LENS  ANTENNAS  MADE  OP  AN  INHOMOGENEOUS 
DIELECTRIC 

Ye.  G.  Zelkin  and  R.  A.  Petrova 

Summary  V. 


,n  ,  v  len^  antennas  ate  examined  with  refraction 

iCtPd  SiJwhcTFT  an?  ^road-directional  lens  antennas  that  are  com¬ 
puted  with  the  help  of  the  conformal  transform  method. 

A  calculation  is  made  of  the  refraction  index  and  the  beam 
trajectory  index  in  these  lenses. 


Introduction 

In  many  practical  cases  it  is  necessary  to  produce  an  antenna 
with  a  broad  beam  pattern,  for  example,  in  a  sector  over  180°.  Develop¬ 
ment  of  antennas  with  these  broad  patterns  presents  great  difficulties. 
It  is  sufficient  to  say  that  the  open  end  of  the  wave  guide  radiates 
only  in  the  120-130°  sector.  In  order  to  expand  the  pattern,  we  have 
to  resort  to  special  measures,  for  example,  the  use  of  dielectric  inserts 
in  the  wave  guide  which  is  constricted  on  the  end.  However,  all  the 
known  methods  of  expanding  the  pattern  generally  impair  the  wide-band 
nature  of  the  antenna. 

The  use  of  lens  antennas  made  of  an  inhomogeneous  dielectric  is 
a  very  promising  trend  in  the  development  of  broad-directional  antennas. 
They  have  a  number  of  advantages.  The  chief  ones  are:  operation  in 
a  broad  band  of  frequencies  and  preservation  of  the  polarization  pro¬ 
perties  of  the* antenna  emitter. 

It  should  be  noted  that  a  broad-directional  pattern  can  be  obtained 
with  the  help  of  lenses  made  of  an  inhomogeneous  dielectric  with  the 
most  diverse  laws  for  the  change  in  the  refraction  index  inside  'f  them. 

This  work  examines  several  lenses  made  of  inhomogeneous  dielectric 
with  a  broad  beam  pattern  (to  360°).  It  calculates  the  refraction 
indices  of  these  lenses,  as  well  as  the  trajectories  of  the  beams 
passing  in  them. 

It  is  common  knowledge  that  the  initial  equation  to  compute  the 
lens  by  the  method  of  geometric  optics  is  the  eikonal  equation 

/  (V  LY=n\  (1) 

where  n--the  function  of  three  variables  in  a  general  case. 

It  is  a  very  complicated  task  to  find  the  solution  to  this 
equation  in  a  general  form  which  satisfies  the  condition  of  obtaining 
a  broad  beam  pattern. 

We  will  examine  below  a  simplified  task  when  the  refraction 

index  in  the  lens  changes  according  to  the  law  'i(A\ «, )  This  law  for 

•\ 

the  change  in  the  refraction  index  makes  it  possible  to  separate  the 
variables  in  equation  (1),  and  to  determine  the  unknown  function  n(<i>). 

The  method  of  conformal  transforms  provides  broad  potentialities 
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for  computing  the  broad-directional  lenses.  There  are  known  lenses 
in  the  literature  that  are  made  of  an  inhomogeneous  dielectric  wh:  ch 
focus  rays  into  parallel  beam  (lenses  of  Mikaelyan ,  Maxwell,  Lyureberg 
et  al.).  The  refraction  indices,  beam  trajectories,  wave  front  and 
other  characteristics  have  been  computed  for  these  lenses.  All  the  rays 
which  emerge  from  the  lenses  indicated  above  are  perpendicular  to  the 
aperture  plane.  '  In  the  case  of  a  broad-directional  lens  ,  the  aperture 
must  not  be  flat ,  but  have  the  shape  of  a  sphere  ,  while  the  rays  must 
emerge  from  it  not  in  a  parallel  but  in  a  spherically  diverging  bundle 
which  is  perpendicular  to  the  aperture  surface.^  Thus,  in  order  to 
obtain  a  broad-directional  pattern,  one  can  conformally  depict  a  lens 
with  flat  aperture  which  focuses  the  rays  into  a  parallel  beam,  into 
a  lens  with  aperture  which  has  the  shape  of  a  sphere.  The  main  pro¬ 
perties  of  the  conformal  depiction:  preservation  of  the  angles  and 
constancy  of  the  strains  will  guarantee  transformation  of  the  flat 
aperture  into  a  spherical  one  ,  and  the  parallel  beam  of  rays  into 
a  spherically  diverging  beam. 

Calculation  of  Lens  Parameters  with  Refraction  Index 

Calculation  of  the  refraction  index.  We  will  examine  a  cylindrical 

lens  with  refraction  index  n  (R,  <p)  -=  liil  where  R  and  <j>--polar  coordi- 

R  • 

nates.  Since  in  cylindrical  lenses,  all  the  beams  lie  in  the  plane 
z=const ,  then  for  simplicity  we  will  examine  the  course  of  the  beams 
in  the  plane  z=0. 

In  our  case,  the  eikonal  equation  in  coordinate  writing  has  the 
following  appearance: 


culation  of  the  lens  we  will  assume  that  it 


has  a  phase 


center  which  coincides  with  the  geometric  center  of  the  sphere.  This 
condition  for  guaranteeing  a  broad  beam  pattern  is  sufficient ,  but  is 
not  necessary  in  the  general  case 
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We  will  look  for  the  solution  in  the  form 


L~  f(R)  -f/.fo).  (3) 

then 

RTW  Jr  /i*(q  )  -  (4) 

Equation  (4)  is  satisfied  on  the  condition 

Rf>(R)~a* 

/;*  (<P)  =  n*  ( Cp)~  flj 

where  a^--random  constant  of  separation  of  the  variables. 

From  here 

f(R)~  f-^dR^OilnR  +  a,  (6) 

where  a--integration  constant: 

fi  (<f)  =  J  Y  nMrf— <P-  '7> 

By  substituting  expressions  (6)  and  (7)  into  (3),  we  have 

L  =  J  j/A2  (<p)  —  a(d<p  -h  Oj  In  R  +  c.  (8) 

l 

One  can  apply  to  equation  (8)  the  study  method  presented  in  boc?k  [1], 
After  differentiating  the  complete  integral  for  a^ ,  and  after  Equating 
the  result  to  the  constant  g^ ,  we  obtain  an  equation  for  the  family 
of  rays  which  depend  on  the  two  constants  a^  and  g^: 

In/?  =  f  -  f1-  —  —  d(f  -f  pi.  (9) 

We  have  obtained  an  integrated  equation  in  relation  to  the  unknown 
function  n(<j>).  The  constant  8^  is  defined  from  the  boundary  condi¬ 
tions,  while  a^  is  the  ray  parameter:  for  each  ray,  the  parameter 
is  a  constant  quantity,  however,  its  value  changes  from  ray  to  ray. 

We  will  determine  the  constant  g^.  For  this  purpose'^  we  will 
find  the  quation  of  rays  which  emerge  from  point  A  (fig.  1),  whose 
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coordinates : 


Figure  1. 


/?,  =  <*,  (10> 


We  will  secure  the  lower  limit  of  the  indefinite  integral  in 
expression  (9)  after  using  the  boundary  conditions  (10): 

In/?  -  f  — ,d +  In d,  (11) 

,aY  «*(<P)-cj 

from  which  it  follows  that  B^=ln  d. 

We  will  now  find  the  limits  for  change  in  a^.  For  this  purpose 
we  will  differentiate  equation  (11)  for 


from  which 


if(lnfi)  _  o i _ 

d<P  \  r  n*  (9)  -  a] 


As  is  known  R-z;  —  lg*, 

dR 

where  a--angle  between  the  radius-vector  and  tangent  to  the  beam  (fig.l) 
We  will  examine  the  rays  emerging  from  point  A  at  angles 


ul<* 


tg “/i=”  corresponds  Co  the  extreme  ray  which  emerges  vertically  up- 
**  n 

wards  ‘P*  ~  2  *  while  tga^=0  corresponds  to  the  central  ray  (<{>^=0). 
Consequently,  the  parameter  a^  which  changes  in  limits  n(Q)l  >  0 
corresponds  to  the  rays  lying  in  the  sector  0  <_ 

We  require  that  all  the  rays  which  emerge  from  point  A  and  which 
are  defined  by  the  values  a^  in  the  indicated  limits  ,  emerge  from  the 
lense  at  points  Bn  which  lie  on  the  circumference  of  radius  R=b.  Then 

K 

f  _2£JL__-  «ln  b.  04) 

fA  }/  d 

Equation  (1A)  can  be  solved  in  relation  to  n(<j>)  if  the  value  n  is 
assigned  at  some  point.  We  will  designate: 

n!(<p)  -  u,  tj>  =  ,  of  =  h.  (15) 

nil 


With  ^q>— 0  u=n2(0)>=«0. 

Since  we  require  that  all  the  rays  emerge  at  points  BR  perpendi¬ 
cularly  to  circumference  R=b ,  then  on  this  circumference  tgcx=0,  and 
consequently,  with 

<P  =  <Pan  ub  =',,(V)  =  0?==/'-  06) 

n  n  n w 


In  the  new  designations  ,  equation  (1A)  looks  like 


n 

f  S~ u 

J 


__  In  b, 

ruT^h  Vh 


where  bQ=  -g,  while  h  changes  in  limits  of  0  <_h<_  Uq.  We  obtained  an 
Abel  type  equation  ,  which  ,  as  is  known  ,has  the  solution 


i  d  l' ,  i  an 

rt  =-  —  --  -  \  In  6fl  -  - . 
*  n  du  J  |  /i  (/:  —  u) 


but  since  ui  =  *** ,  ,  then 
du 


(18) 


where  A 


After  taking  the  integral  in  the  right  side,  we  obtain 

*  “  A  ln  “7  V— - •  (20) 

2  u®  -  j.  2 u#  —  u 

By  solving  (2)  in  relation  to  u,  after  certain  transforms  we  obtain 

u  — - — — — —  .  (21) 

«) 

2  2 

but  since  u=n  (<j>)  ,  while  Ug=n  (0)  ,  then  the  final  expression  for  the 
refraction  index 

(22) 


We  see  that  the  refraction  index  on  the  lens  aperture  depends 
on  the  selection  of  n(0)  ,  as  well  as  on  the  arrangement  of  the  focus 
in  the  lens.  It  is  assigned  by  the  quantity  d.  The  dependence  of 
the  refraction  index  on  the  lens  boundary,  i.e.  ,  with  R=b  ,  on  the 
selected  parameters  n(0)  and  d  is  presented  in  the  table  below.  It 
does  not  present  the  absolute  values  of  the  refraction  indices  ,  but 
the  ratios  of  the  refraction  indices  with  different  <J>  to  n(0).  The 
calculations  were  r : ie  for  d=0.1  and  0.2b. 

It  is  apparent  from  the  table  that  the  greatest  value  of  the 
refraction  index  is  with  4>  =0 .  With  an  increase  in  <j>  ,  the  refraction 
index  diminishes ,  as  the  change  is  larger  the  greater  the  quantity  d, 
i.e. ,  the  further  the  focus  is  located  from  the  center  of  the  circle 
which  is  the  lens  aperture. 

TABLE  _ _ _ 


d  in  relative 
units 


n(<f>) 


Values  of  relative  quantity  — with  different  <|> 


30 

60-  j 

; 

120’  | 

I50: 

ISO’ 

0,94 

0,79  1 

0,61 

1 

0,43 

0,33 

0,23 

0,88 

0.64  ] 

i 

0,41 

0,26 

0,15 

0.1 

The  complete  refraction  index  in  the  lens  can  be  expressed  by 
taking  into  consideration  (22): 


« (V. 


R)  - 


Sch 

..  b 


(23)  ^ 


Calculation  of  the  ray  traj  ectories .  In  order  to  find  the  ray 
trajectories,  it  is  necessary  to  substitute  expression  (22)  into 
formula  (11)  and  to  compute  the  obtained  integral: 


<p 


yr_2!m__a* 

V  K)' 


(24) 


We  will  replace  the  variables:  x=shk<{> ,  k  =  -*-r  and  cs=n,'°>( 
then,  equation  (24)  adopts  the  appearance 

Af  »h  x 


lnf=  f  —±r . - 

d  ,)  V  et -*»- 


(25) 


The  right  side  of  expression  (25)  is  the  tabular  integral  [2]. 
By  solving  it ,  we  obtain 


j  1  • 

In  -•  =  —  arc  sin 
d  k 

Arshx 

A/iTZT  ’ 

(26) 

from  here , 

.  R 

„  T  '  «(0)]*  , 

*  V  m-1 

- .  to  T 

Sm  2 A  ' 

(27) 

Formula  (27)  determines  the  link  between  $  and  R  with  the  selected 
values  of  d  and  n(0).  By  assigning  the  parameter  a^  in  limits  from 
n(0)  to  0  ,  we  will  obtain  the  trajectories  of  the  rays  which  emerge 
from  the  focus  at  angles  from  0  to  90°. 

The  trajecf~-*"\  .  fo”  e  rays  for  the  two  cases  where  d=0.1  and 
0.3b  are  presented  in  fig.  2  which  only  depicts  the  upper  half  of  the 
lens.  The  entire  lens  in  this  case  has  the  appearance  of  a  sphere 
with  removed  middle,  where  the  emitter  can  be  placed. 

With  d=0.1b,  we  are  c  _  limited  to  constructing  the  rays  that 
emerge  from  the  focus  at  angles  from  ^=0°  to  74°  [in  this  case,  the 
parameter  a^  changed  from  n(0)  to  0.23  n(0)].  The  outer  ray  (4>i=74°) 
emerges  from  the  lens  perpendicularly  to  the  aperture  at  angle  4>  =1 30°. 
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With  d=0.2b,  a  ray  emerges  from  the  lens  at  angle  <f>=180°.  It 
travels  from  the  focus  at  angle  ^=71°.  This  corresponds  to  the  value 
a^  which  equals  0.1  n(0). 

Depending  on  the  selection  of  values  n(0)  ,  one  can  obtain  beam 
patterns  of  varying  width.  With  d=0.1b,  for  example,  the  refraction 
index  on  the  lens  aperture  changes  from  n(0)  to  0.23  n(0)  with  <j>  =tt. 

This  means,  that  with  n(0)=4.3,  we  can  create  radiation  in  the  360° 
sector.  With  a  value  of  n(0)=1.6,  with  the  same  d=0.1b,  the  emitter 
ray  emerges  from  the  focus  at  an  angle  of  <^=52°. 

Lens  Calculation  by  the  Method  of  Conformal  Transform 

Selection  of  the  parameters  of  the  transforming  function .  We 
will  first  examine  the  known  lenses  of  Maxwell  and  Mikaelyan  [1]  as 
the  initial  lenses.  In  the  Maxwell  lens,  the  refraction  index  depends 
only  on  the  radius  of  the  lens  ,  while  ir.  the  Mikaelyan  lens  it 
depends  on  one  coordinate  y.  If  we  examine  both  lenses  in  the  Cartesian 
coordinate  system,  then  the  aperture  of  both  lenses  will  be  a  plane 
perpendicular  to  the  x-axis  and  intersecting  it  at  point  x=0. 

We  will  examine  the  pattern  of  rays  in  one  plane,  for  example, 

2=0.  All  the  rays  in  the  indicated  lenses  intersect  the  y-axis  at  a 
right  angle.  In  order  to  obtain  a  broad-directional  lens,  the  rays 
must  intersect  the  circle  at  a  right  angle ,  consequently ,  this  lens 
can  be  produced  if  the  y-axis  is  conformally  transformed  into  a  circle. 

We  will  use  the  linear-fractional  function  as  the  transforming 
function 


H3 


where  U"~u  +  ip,  ? -rx+\y.  ,  while  a,p,  x  and  y--resocctively  the  coordi¬ 
nates  of  the  broad-directional  and  transformable  lens.  For  simplicity, 


Figure  3. 

one  of  the  coefficients,  for  example,  k,  can  be  assumed  to  be  equal  to 

a  unit  [3],  Thus,  the  illustration  is  assigned  by  three  coefficients 

1,  m  and  n  which  depend  on  three  complex  or  six  actual  parameters. 

< 

We  will  first  examine  the  Milcaelyan  lens.  Assume  that  the  lens 
aperture  will  be  2b  (fig.  3a).  We  will  require  that  the  refraction 
index  equal  one  along  the  straight  lines  y=±b.  We  will  transform 
this  aperture  into  part  of  a  circle  which  is  defined  by  the  coefficients 
b(a^±iB^)  (fig.  3b).  Since  it  is  sufficient  to  assign  three  points 
in  order  to  define  the  circle  ,  we  will  take  the  following  as  the  original 
three  points  on  the  lens  aperture:  Z=ib,  0,  -ib- ,  and  we  will  transform 
them  into  the  points  W=b(a^+ig^)  ,  b,  b(ct^-iG^).  We  will  assume  for 
simplicity  that  the  quantity  b=l.  In  order  to  find  the  coefficients  1, 
ma  and  n,  we  substitute  the  values  W  and  the  values  of  Z  corresponding 
to  them  into  formula  (28): 


n 


i  \l' 
i7i  i  +  n 


*i—  ipi 


—  i  4-  < 
—  m  i  +  n 


(29) 
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From  which,  after  certain  transforms,  we  obtain  the  values  for  all  three 
coefficients 

n  ~  /  —  — ■ —  - 

(3C) 

m  «  —  1 


After  substituting  the  obtained  values  for  the  coefficients  into 
formula  (28),  and  after  separating  the  real  and  imaginary  parts,  we 
obtain  an  expression  for  x  and  y  through* the  coordinates  a  and  g: 

/l(a»t  Pl_-  h_  (31) 

Kl  +  aP+W 

y  = - ML—  ,  (32) 

y  1(1  +  a)*+  P'J 

where  A  —  — — -  . 

1  — <*i 

It  is  easy  to  see  from  formula  (31)  that  the  y-axis  ,  which  the 
rays  in  the  Mikaelyan  lens  are  perpendicular  to  ,  becomes  a  circle 
which  is  described  by  the  equation  a2+g2=  1.  This  was  also  required 
in  the  statement  of  the  problem. 

Refraction  index  and  trajectories  of  rays  in  a.  broad- directional 
lens  of  the  Mikaelyan  lens  type. .  It  is  common  knowledge  II]  that  in 
a  Mikaelyan  lens  of  single  radius  ,  the  refraction  coefficient  only 
depends  on  one  coordinate  y: 


n  (</) « 


In  order  to  guarantee  invariance  in  the  optical  path,  the  refraction 
index  in  the  new  lens  must  be  transformed  according  to  the  formula: 

w(»,frt=-=i‘§7f-.  l*> 

where  n(a  ,8)  can  be  obtained  from  formula  (33)  after  replacement  of 
y  through  the  new  coordinates  according  to  (32),  and  jW'  | 
modulus  of  derivative  from  analytical  function  W(z).  By  using 
expressions  (28)  ,  (30)  ,  (31)  and  (32)  ,  we  find 

<»> 
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Consequently,  the  final  expression  for  the  refraction  index  in  the 
broad-directional  lens  for  a  Mikaelyan  type  lens: 


A'  (*,  p) 


2|,1  !.•)<<') 


|(1  4-  a)»  i  P’Jch 


:i.  lP  __ 

(1  :  u>«  :  P* 


(31a) 


Formula  (34a)  makes  it  possible  to  compute  the  refraction  indices  in 
a  broad-directional  lens. 

The  ray  trajectories  in  the  Mikaelyan  lens  of  a  single  radius  are 
described  by  the  equation 


. . _ «i . s!l  •*>  ft 

2  \'r  «5  ('!)  —  2 


(3t») 


where  a^--paratneter  which  is  constant  for  each  ray  and  changes  from  ray 
to  ray.  For  any  ray  in  a  bundle  0  <  (|ij<  ^  ,  the  parameter  a^  satisfies 
the  condition  0  ai  1  n(0). 

After  substituting,  instead  of  x  and  y,  expressions  (31)  and  (32), 
we  obtain  an  expression  which  determines  the  ray  trajectories  in  a 

broad- directional  lens: 

\  . 


cos  ii?  — 

2  ((l+ap-pq 


—  — 0'_ _ — .  sh  — 

J  '  1(1 


n  .1 P _ 

r«P  p=j  ■ 


(37) 


The  focus  of  the  Mikaelyan  lens  z^  which  has  coordinates  x^®  -1, 
y^,s0 ,  according  to  the  formulas  for  transformation  of  the  coordinates , 
becomes  point  w.  with  coordinates: 


Pi  +  a*  ~  1 
1  —  <>i  +  Pi 


P*  -■  0. 


(38) 


Formula  (38)  shows  that  the  position  of  the  focus  in  the  transformed 
lens  is  only  determined  by  tha  position  of  the  extreme  rays  in  the 
lens.  In  a  particular  case,  where  the  sector  of  the  rays  equals  ±90°, 
the  lens  focus  coincides  with  the  geometric  center  of  the  circle. 

A  lens  with  sector  ±120°  was  examined  as  an  example.  For  it, 
a^  =  -0.5,  $2=0.866,  A=  -0,577.  In  this  lens,  the  focus  is  located 
at  the  point  with  coordinates  B$=0,  oi^®  -0.27.  The  lines  of  equal 
refraction  indices  (dotted  line)  and  ray  trajectories  in  the  lens 

(solid)  are  depicted  ir.  fig.  4. 
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W:e  examined  above  the  transform  of  a  Mikaelyan  lens  of  single 
radius  in  a  broad-directional  ltns,  also  of  single  radius,  and  it  was 
hypothesized  that  the  refraction  coefficient  equalled  1  with  y*±l.  One 
can  examine  a  more  general  case  where  the  lens  size  equals  b,  while 
n(yW  with  a  certain  y=a.  The  aperture  of  this  lens  equals  2a,  and  the 
quantity  a  is  associated  with  the  refraction  coefficient  n(0)  by  the 
correlation 

a  a  -^archn(O).  (39) 


The  coordinate  transform  formulas  in  this  case  look  like 


aA  (a*  +  p_-  >*)_  20  M _ 

X  —  :  v  «...  1  »  ((a  -{•  b)*  +  PM 


The  refraction  coefficients 


(40) 


#»(*.««  - 


n(Q) 


,  ,  aM 
r,'d>  l«  +  *)«+M 


(41) 


The  trajectories  of  the  rays  are  described  by  the  equation 


cos  :i aA  _j£  hPinfc!_  _ ft _ _  sh  nfl4p _ 

2  b  (o  +  6)«  +  P»  (o  +  6)»+P> 


(42) 


Figure  4. 

Refraction  Coefficient  and  Traj ectories  of  Rays  in  Broad- 
Directional  Lens  of  the  Maxwell  Lens  Type .  We  will  now  examine  the 
Maxwell  lens,  for  example,  of  single  radius  and  will  transform  it 
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into  a  broad-dircctional  lens  also  of  a  single  radius.  This  transform 
is  reduced  to  the  previously  examined  transform  of  the  Mikaelyan  lens 
with  aperture  equal  to  2  ,  and  the  refraction  coefficient  equal  to  1 
with  y=±l  ,  and  consequently,  the  transform  coefficients,  coordinates 
x  and  y,  and  the  modulus  of  the  derivative  W(z)  are  reflected  corres¬ 
pondingly  by  formulas  (30),  (31),  (32)  and  (35). 

In  a  Maxwell  lens  of  single  radius,  the  refraction  coefficient 


n(x,  </)« 


MO) 

i  4“  Jt*  -f  vs 


(43) 


while  in  the  broad-directional,  obtained  on  the  basis  of  this  lens, 
correspondingly 


A'(>,  P) 


_ 2 1  A  !  n  (0) _ 

Wl  -i-  aj4  P4)  (I  -|-  ,1=)  -  ia  ‘ 


(44) 


It  is  easy  to  show  that  the  lines  of  equal  values  for  the 
refraction  coefficients  N(a,e)=const  are  circles  with  center  that  is 
shifted  in  relation  to  the  beginning  of  the  coordinates.  These  circles 
are  determined  by  the  expression 


where 


6=  1  +  a,  a 


r_  N(a,  P) 
MO)  ' 


2(2>1»C— Ml) 
+ 1)  ’ 


0, 


(45) 


The  center  of  these  circles  is  located  at  the  point  with 
coordinates : 


a 


o 


A*+  1  ’ 


Po-0. 


(16) 


while  the  radii  of  the  circles 
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The  lens  focus  is  located,  as  before,  at  the  point  with  coordi- 


Pt-f  0|  — > 


0. 


nates : 


a. 


The  trajectories  of  the  rays  in  the  Maxwell  lens  of  single 
radius  are  described  by  the  equation 

x*  -b  (/*  +  2«/etgy  -  1,  (  IS) 

where  Y--angle  of  emergence  of  the  ray  from  the  lens  focus. 

After  substituting  expressions  (31)  and  (32)  into  formula  (48), 
we  obtain  a  formula  to  compute  the  trajectories  of  the  rays  in  a  broad- 
directional  lens: 


/}»((,» 4.  pa  -  1)5. t.  4A  P  clg  y 

i(i  +«)*  +  P:]3  '  id +®is~i-P2i 


(19) 


After  transforming  this  formula  ,  it  is  easy  to  be  convinced  that 
the  trajectories  of  the  rays  have  the  shape  of  circles  which  depend 
on  the  y  parameter: 


A*  -H 
A-—  l 


2»  .Mdfi.lf.  +  l 

1  A*  -  1 


0. 


(19a) 


with  centers  at  the  points  with  coordinates: 

„  _  A-  -|-  1  o  2.4  dg  v 

0  A'--  1  ’  ™  A* -l 


(50) 


and  radii 

As  an  example,  the  lens  with  radiation  sector  ±120°  was  computed. 
The  lines  of  equal  refraction  coefficients  (dotted  lines)  and  the 
trajectories  of  the  rays  in  this  lens  (solid  line)  are  presented  in 
fig.  5. 

Calculation  of  a  broad- directional  lens  according  to  the  type 
of  isotropic  emitter .  The  method  of  conformal  transformation  makes  it 
possible  to  obtain  broad-directional  lenses  not  only  based  on  the 
known  lenses,  but  also  quite  new  ones. 

As  an  example  ,  we  will  examine  a  broad-directional  lens  which 
was  obtained  based  on  an  isotropic  emitter  illustrated  in  fig.  6. 

It  is  assumed  here  that  the  hypothetical  isotropic  emitter  (fig.  6a) 
is  located  in.  the  center  of  a  spherical  lens  made  of  a  uniform 
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Figure  5. 


dielectric  or  nonuniform  dielectric ,  but  with  refraction  index  that 
changes  only  for  the  radius. 

As  before  ,  we  will  examine  the  pattern  of  distribution  for  the 
refraction  coefficients  and  the  ray  trajectories  in  the  Z=0  plane. 

We  will  conformally  depict  the  right  semicircle  with  outer  rays  at 
points  with  coordinates  x=0  and  y=±l  into  a  broader  spectrum  with 
outer  rays  at  the  points  with  coordinates  (fig-  6b).  The 

unknown  piecewise-linear  function  for  this  transform  has  the  coeffi¬ 
cients  : 

,“m="r+E  m: 

t:  ' 

and  is  written  as  follows: 


where 


l  +  Pi ' 


£*  + 1  * 


<53) 


By  separating  the  actual  and  the  imaginary  parts ,  we  obtain 

expressions  for  x  and  y  through  the  new  coordinates  a  and  6: 

*  q(l  -t-  /!•)  —  a fn*  ■  ft-  i) 

(I  —  u/J)-  -i  /i* 

_ ft(i  -/i) 

(I  -  ally1  -j  y*li‘ 
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(65) 
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Figure  6. 


If  we  assume  that  the  refraction  coefficient  in  the  original  lens 
has  a  constant  value  n(0),  then  we  can  find  the  refraction  coefficient 
in  a  broad  directional  lens  ,  after  preliminarily  determining  the 
modulus  of  the  derivative  W(z) ,  For  this  transform 


1 W I  - - 


1 1  -  fl*  l 

</*+  \r--viihji 


where  x  and  y  are  determined  by  formulas  (54)  and  (55).  The  final 
expression  for  the  refraction  coefficient 


AM*.  P) 


w (0) (i  _ b) in  -l  ann-a^)3 

(1  +  #)  I(l  '—o  W  P: 


(57) 


The  point  of  location  of  the  emitter  whose  coordinates  are  x=0  ,  y-0 

in  the  transformed  lens  becomes  the  point  with  coordinates 

B#-0. 

In  the  original  lens,  the  ray  trajectories  are  segments  of 
straight  lines: 


lev. 


(58) 


where  y  is  the  angle  of  the  ray  emergence. 

In  the  new  lens,  the  ray  trajectories  are  circles  described  by 
the  equation 
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modeling  beam  patterns  of  linear  antennas 
K.  s.  Shcheglov 


Summary 


A  method  is  examined  for  mathematical  modeling  on  a  computer 
of  the  beam  patterns  of  antennas  based  on  the  use  of  Kotel 'nikov’ s 
theorem.  The  values  of  the  patterns  for  a  number  of  reference 
directions  are  used  for  complete  analysis  of  the  properties  of  the 
pattern,  including  directive  gain  and  scattering  coefficient. 

The  examined  algorithm  is  applied  to  computing  antennas  by 
the  method  of  statistical  tests  (Monte  Carlo  method)  which  con¬ 
sists  of  modeling  random  patterns. 


Introduction 


Difficulties  in  analytical  examination  of  beam  patterns  of 
antennas  with  complicated  amplitude-phase  distributions  in  the  aperture 
have  led  to  the  need  for  development  of  special  methods  of  study.  They 
include,  on  the  one  hand,  the  methods  of  graphic  modeling  of  the 
patterns  using  different  types  of  analog  devices  and  plans  equivalent 
to  the  antenna  [1].  On  the  other  hand,  the  development  of  computer 
equipment  has  afforded  broad  possibilities  of  using  the  methods  of 
mathematical  modeling  of  antenna  patterns  [2,3].  As  compared  to  the 
methods  of  the  first  group  ,  these  methods  reduce  to  a  minimum  the 
loss  of  time  and  resources  with  considerably  greater  accuracy  of  the 
results. 

Pattern  modeling  is  especially  important  in  studying  antennas 
with  random  distributions  in  the  aperture  caused  by  errors  in  fabri¬ 
cation  or  which  are  especially  assigned.  The  statistical  theory  of 
these  antennnas  [4]  makes  it  possible  to  define  the  characteristics 
for  the  main  pattern  parameters  only  in  a  statistical  sense.  Only 
the  class  of  normal  errors  is  subject  to  detailed  analysis.  These 
characteristics  determine  the  properties  of  the  pattern  in  a  group 
of  single-type  antennas. 

It  is  important  to  directly  calculate  the  individual  realiza¬ 
tions  of  random  patterns  and  to  make  a  subsequent  statistical  study 
of  the  parameters  of  the  pattern  group.  These  tasks  are  solved  in 
the  framework  of  the  method  of  statistical  tests,  the  Monte  Carlo 
method  [5]  which  has  become  popular  with  the  development  of  computer 
equipment.  In  the  most  general  form,  the  Monte  Carlo  method  consists 
of  solving  different  problems  of  computer  mathematics  by  constructing 
a  random  process  for  each  case  with  parameters  which  are  equal  to  the 
unknown  quantities  of  this  problem.  As  applied  to  computing  antennas  , 
this  means  that  a  number  of  random  realizations  are  obtained  on  the 
computer  for  distribution  in  the  aperture  with  required  statistical 
properties  and  the  patterns  corresponding  to  them  are  defined. 

Based  on  a  number  of  realizations  of  random  patterns,  the  proba¬ 
bility  estimates  of  their  parameters  are  determined,  including  those 
such  as  coefficients  of  scattering  which  are  difficult  to  obtain  by 
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other  methods.  In  this  case  the  law  of  distribution  of  errors  can  be 
arbitrary. 

The  method  examined  below  can  be  applied  to  problems  of  heuristic 
calculation  of  antennas.  In  this  case,  the  modeling  method  is  used 
to  select  the  most  successful  random  distributions  of  discrete  emitters 
(or  their  amplitudes)  which  yield  satisfactory  patterns.  A  further 
task  in  this  case  is  accurate  fulfillment  of  the  found  random  function 
for  the  arrangement  of  the  elements  (distribution  of  amplitudes). 

We  note  in  conclusion  that  in  addition  to  the  tasks  of  computing 
antennas,  modeling  the  patterns  is  also  used  to  model  systems  in  which 
these  antennas  are  used. 


Discrete  Presentation  of  the  Beam  Patterns 

Assume  that  in  the  aperture  of  a  linear  antenna  of  length  L, 
the  distribution  a(s)  is  assigned  and  is  equal  in  a  general  case  to: 

'  a  (s)  =  a0(s)  +  y,(s)  -j-  i  *(s).  (1) 


The  beam  pattern  of  this  antenn 


f(u)  =  jaG)e"*<n, 


(2) 


K  l  A 

where  u=  ^-cos  t>,  V~ --angle  which  is  read  off  from  the  antennna  axis, 
5--normed  coordinate  of  the  aperture  point,  k--wave  number  in  a  free 
space  (fig.  1). 

If  the  antenna  aperture  is  stimulated  by  a  wave  with  wave  number 
Y,  one  should  place  u-8  in  expression  (2)  instead  of  the  quantity  u 
in  the  exponent,  where  B  =y  This  replacement  does  not  affect  the 
nature  of  the  F(u)  relationship  and  must  only  be  taken  into  considera¬ 
tion  in  the  transition  to  the  angular  coordinate  in  particular  ,  in 
computing  the  directive  gain  and  the  scattering  coefficient  of  the 
antenna. 

In  accordance  with  (1)  ,  we  present  the  pattern  (2)  in  the  form 


F  (u)  =  F0  (u)  +  Fr  (w)  -I-  i  Ft  («).  (3) 
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Figure  1 . 


If  in  the  distribution  of  a(s)  we  can  isolate  a  component  which 
is  convenient  for  analytical  integration,  it  should  be  included  in  aQ  , 
otherwise  ,  we  consider  3-q=0  and  correspondingly,  Fa=0. 

For  convenience  we  will  introduce  the  following  designations: 


;  i 

It  —  l  Ms) sin d  |  /(g) sin «  Idl.  (4) 

-»  1 1 

Thus,  the  pattern  (3)  can  be  presented  in  the  form 
+  (£“<*)  -  i(fc  -1-  It).  (5) 

Direct  numerical  computation  of  the  pattern  with  the  help  of 
formula  (2)  is  inconvenient  with  a  large  number  of  points  for  u.  In 
addition,  determination  of  a  number  of  important  characteristics  of 
the  antenna  (amplification,  scattering  coefficient,  position  of  the 
main  maximum,  etc.)  by  direct  methods  is  associated  with  cumbersome 
calculations . 

In  this  respect,  we  will  use  presentation  of  the  pattern  in  the 
form  of  a  series  for  reference  functions  that,  for  etample ,  are  very 
popular  in  the  tasks  of  antenna  synthesis  [6,7].  According  to  Kotel'- 
nikov's  theorem  [7,8]  or  the  theorem  of  readings,  the  beam  pattern  examined 
as  a  function  with  limited  spectrum  (distribution)  can  be  presented 
in  the  form  (fig.  2)  of  a  series  for  the  functions  S1-^-  —  : 


F(u)  --  \ i'fta,,)^) 
(u  --  m  n) 
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The  reference  functions  - —  are  arranged  uniformly  with 

interval  Au=?.  on  the  u  axis.  The  coefficients  in  (6)  with  them  are 
equal  to  the  values  of  the  pattern  F(u)  in  the  reference  points  u- 
nuu  ,  m=0  ,  ±1  ,  ±2 .  .  . 

"We  will  isolate  the  real  and  imaginary  parts  in  pattern  (5): 

F  =X-j-\Y, 

•where  X  -  j40  (m)  ->•  Ai  («);  Y  -■  C0 («)  +  Ct  ( u ), 


s*nd  Aj(u)  and  C^(u)  according  to  (6)  equal: 

A  <«)  =  V  , 

jimA  (u  —  m  n) 

c,c O  ' 

jm J  (ii-bij) 


(8) 


If  we  con  Isolate  the  component  aQ  in  distribution  (1)  ,  the 
function 

i 

A>(m)  -  j’  a0(l)cosuld  l  (9) 

I  .correspond  to  the  "summary"  pattern,  while 

i 

C»  (")  /  °o(5)sinw“dt 

ho  "differential . " 
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The  index  m  in  the  functions  8m"clm  in  (8)  means  that  they  were 
computed  for  the  value  of  the  argument  u=mTt .  The  functions  g  and  b 
are  even  for  u;  h  and  d  are  odd.  Expressions  (8)  can  therefore  be 
presented  in  the  form: 


A  (±  «)  =3  ^  'm  1  Sm  (sm  +  S— ml  ~  (S"  S~m) I  ’ 

m— 0 

%m  (Sm  +  s-«)  ~  hm ~ s~m^ * 


-i-  withm  -  0, 

2 

lwith  m  —  1,2,  . 


The  following  designation  is  adopted  in  formulas  (11) -(12)  for 
convenience 

sin  (ti  ~~m  n)  /10v 

»m - —  .  I  *  o; 

v«  —  m  n) 

In  the  reference  directions  u=mn ,  the  function  sm~l ,  the  remaining 
sn=0  (n#n)  and  the  function  F(u)  in  (6)  equals  its  reference  value. 

Directive  Gain  and  Scattering  Coefficient 

By  using  the  discrete  concept  (6)  of  the  beam  pattern,  we  can 
find  expressions  for  the  directive  gain  which  are  convenient  for 
calculation  on  a  computer 

c  r.  if  e-iKt.  04) 

MS 


and  scattering  coefficient 

"  c  -  - 1 


of  a  linear  antenna. 

In  expressions  (14)  and  (15)  ,  Pma^c  is  the  level  of  power  in  the 
main  direction,  P  --complete  emitted  power,  and  PrAJl  --powers 
emitted  in  the  region  of  lateral  lobes  and  in  the  main  lobe  respectively 
For  a  linear  antenna  ,  the  complete  emitted  power 

/  S'  *  «»  ... 

' .  Pr  **  f  (  p («) sin •»  d !> d «r  I /»(«)</«,  {\G} 


Figure  3. 


kL  kL 

where  u^*= - (1  +  A)  and  U2~  — j  (l-A)--boundaries  of  change  in  the 

generalized^variable  u  which  correspond  to  the  boundaries  of  the  real 
pattern  lV=0  and  it,  while  A  =  -y/k. 

The  power  which  is  contained  in  the  main  lobe ,  is  determined 

by  the  same  correlation  (16)  with  replacement  of  u^  and  U2  by  the 
coordinates  u-^q  and  U2Q  (fig.  3)  of  the  boundaries  of  the  main  lobe. 

If  A  >  0,  the  main  lobe  is  inclined  towards  the  axis,  therefore 
its  right  boundary  i^a^c  is  determined  from  the  condition 


Wlo  ^ 

^  ]Q  Wj, 


Since  the  expression  for  the  pattern  p(u)  consists  of  three  terms 


P  («)  -  (Al  +  C|)  f  2  (-Vh  4-  C0C,)  4-  (A*  +  CJ).  (18) 
one  can  correspondingly  divide  the  integral  (16)  into  parts 

nx  --  I'  p(u)du  -  nx+  /7,+  77*.  (19) 

We  will  first  compute  the  last  terra  in  the  right  side  of  (19) : 

/V-  ]  (A\+C])dur 
— «» 

K  V  v  10JW  -  dm)  (ft.  -  d„)  +  (pm  +  K)  (>>n  +  K)\  Tm,  „  •  (20) 


in  w  n— — •> 

The  coefficients  in  series  (20) 


T  (~Ul  u,)  -  f  ^^LZ^Ldu. 

m. " '  1  J  (u  —  m  n)  (u  —  n: 


(21) 
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In  the  infinite  interval  «*  <  u  <  ,  the  functions  s  [see  (11)] 

m 


are  orthogonal  with  m  whole 

7V  n  ’  where 


l  m  «  n 
0  m  ft. 


For  finite  values  u,  and  u, ,  the  functions  T_  _  equal 

A  4  IU  |IT  % 


rr>. «,)  -  z.  (o)  ;  ~L( v)  j  .  (22) 

'  — *l— -Hfi— UK  / 


The  functions 


L  (v)  — -  In  v  —  Ci  (2u) 


are  even.  This  follows  from  presentation  of  the  integral  cosine  Ci(x) 
indie  form  of  a  series,  therefore  the  limits  of  integration  in  (22) 
are  taken  according  to  the  modulus 


T*.  «<-«».  “*)“Q(p) 


•where 


<?(t)  -  Si(2y)  — 


1  —  cos  2d 


The  function  Q(v)  is  odd. 

For  the  functions  L(v)  and  Q(v)  ,  one  can  obtain  expressions  in 
form  of  the  series : 

L(v)  —  V  (—  1  )*~*  -u2~- ;  Qlv)~  1  )*~1  - (24) 

JJT  m  (2k)  \  ’  (2A-  —  1)  (2Ji)!  ^  1 


the  function  L(v)  ,  the  constant  c-ln2  which  is  not  important  for 
<£^->p.utation  of  T  is  omitted,  where  c--Euler' s' constant)  .  With 
v  •>>  J,  one  should  substitute  in  L(v)  and  Q(v)  the  asymptotic 
^prjessions  for  the  functions  Si(x)  and  Ci(x). 

After  using  the  properties  of  symmetry  of  the  functions  (A)  for 
end  the  properties  of  the  functions  n=Tn  m»  one  can  significantly 
^\eduee  the  double  sum  (20)  which  is  cumbersome  to  compute,  by  pas&ing 
iV  summation  in  the  sector  0<_  m  £  ,  m  <^n  <  «. 

The  final  expression  for  II^:  * 
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n»  =  V  ^  Fm>  H  *„*•  |(M„  :  Until)  Ti  4 


,  (25) 


4  (b  A  —  Km<(„)T,  4  (bnhm  —  y„Jm)  Tj  4  (/«,A  4  </„<(„)  T4 


where 


lwithm^n,  #  =  1 

2w<thm^/i;  m  1 


1/2  withm  - 
*  with  M  T“  (.  2,3  ,  .  . 


The  functions  t  equal: 


'  ti  *=  T  n-hTm^,  4  T_n  .  4  T_„  •  t»  -•  Tm  „  4 

*  ffij  n  #Fii*— #1  — /n*  ^  1  n*  nt%  n 

4-  T  —  T  —  T 
r  m»— n  #  — /M,n  — m.—** 

.»  --  T  —  T*  4-  T*  _ T  •  r  T*  _ .  T  _ 

*•  *  m,  n  1  m,— n  ‘  1  -hi,  n  1  «•  *4  1  m,  n  9  m,— it 

-r.„.„+r — ,  (20) 

With  discrete  presentation  of  the  patterns  (6)  ,  there  is  no  diffi¬ 
culty  in  computing  the  first  two  terms  in  (19).  For  this  purpose  we 

can  use  direct  integration  for  u.  If  the  pattern  F0=A04iC0  can  be 

"  s  in  x 

represented  as  a  finite  number  of  terms  of  the  type  of  functions  — - — , 
the  integral  n<—2  j  (4oAt4C(4':)</ii  is  expressed  by  a  simple  series. 

Thus  ,  if  in  the  distribution  we  isolate  the  component  of 
"summary"  distribution  of  the  type  ao8*=l+4n*'  ,  the  pattern 

looks  like: 

A,  *■  so  +  y  7  (si  +  s_t);  cK  -  o. 

In  this  case 


"»  -  \  o  +  9 (**..,  H-  T'o.  -.)  +  V  W. «+  27V ->  +  r->.  -« ) 

m 

ft*  --  2  |£m  4  P_m)  dm  (j'm  !’m)J 


,  (27) 


where 


Vm  J  A°- Smdu  ro.  m  +  y  <7  (7’It  m  4  T_,  m). 


61 


For  the  "differential"  distribution  ~  c^sin  psin  ”  ^ 

the  pattern 

=*  0;  Ctl  ----  c |  (s,  *  ip/Sj-s  ,M 
{  Vs  *J J 


^  ^  j(r..  |-~7'i.-i)+  2p^r  L-  r  _L^  + 

+  PV7J  j-7«  _A)  *  , 

v  2  2  2‘  v)  i 

m 

H*  =*  2C  ^  £*  !  *">  (V*  +  *_m)  +  Am  (vm  “  v_m)J  I 


vm  *  J  C0i smdu  —  (TUm  m)  4-  p  I^T j  ^  —  T_  j  ( m  V 


The  functions  Tm  r  with  semiwhole  indices  equal: 

/  u,-m-  */2 


— s/2 


r  !C_  SK20)  |  +SUS0) 

"•  *  J-  I  m  T  ~^j2a  l 


_2_ 

J._J  n 
2  2 


=  —  L(2i>)  —L(2v) 


~u‘+l 


% 

~~Ut  2 


It  is  easy  to  find  similar  expressions  in  other  cases ,  if  the 
amplitude  distribution  of  aQ  can  be  expanded  into  sinusoidal  components 
If  this  is  difficult  to  do,  Sq  may  not  be  isolated  in  distribution  (1). 
Then  and  n 2  automatically  equal  zero,  since  the  pattern 

^0  —  A>  -f  i  C,. 


equals  zero. 


Thus,  we  finally  obtain  expressions  for  the  directive  gain 


G~kL- 


Pui kc (u) 


IFor  calculation  of  G  we  can  use  a  very  convenient  method  which  is 
presented  in  [8],  however  it  is  not  suitable  for  computation  of  t ■ 
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^**.*aj  ^  1 ii—  ^  j.iw ~i i>  •  %«>.  u 


S' 


<>■  •*'' 

i> 

i! 


r  * 


and  the  scattering  coefficient 

n,, 

fly  ' 


t-  1 


(30) 


The  power  which  is  contained  in  the  main  lobe  ii^q  is  determined  by 
the  same  ratios  as  n ^  ,  with  replacement  of  u-^  and  by  the  coordinates 


of  the  first  minimums  of  the  pattern  u^q  and  U2q 


[see  (17)1. 


Modeling  the  Patterns 

The  effectiveness  of  mathematical  modeling  significantly  depends 
on  how  economically  the  caclulation  program  is  compiled.  In  the 
examined  case  of  the  discrete  presentation  of  the  pattern,  it  is 
important  to  select  M,  the  number  of  terms  in  Kotei ‘ nikov 1 s  series 
(11)- (12)  which  yields  a  sufficient  degree  of  approximation.  As  shown 

90 

in  [7]  ,  the  standard  error  of  approximating  the  function  F(u)  ~  V  F(k  n)  sk 

by  the  finite  series  (u)  —  ^ ^  (k  a)  sk  equals  the  energy  of  the 

rejected  "tails"  in  the  series  n{  V  |f(A.i)  |* -f 


Af-H 


If  the  F(u)  functions  diminishes  as  — ^,then  one  can  find  an 

v 

estimate  for  the  approximation  error: 

l  F  («)  ■ -  («)  S  < - “Vr  •  W) 


a*-] 


At  the  points  u=nw  ,  the  error  <$m  equals  zero. 


From  condition 

(31)  ,  we  compute  the  number  M  which  yields  permissible  error  at 
points  u  =  m-rr  +  %  which  also  replaces  the  infinite  limit  of  summation 
in  the  formulas  of  the  preceding  sections. 

Modeling  of  the  pattern  begins  with  computations  for  the  selected 
amplitude  distribution  (1)  of  the  reference  (u=m^)  values  of  the 
functions  Sm-dm  for  formulas  (A)  for  0<m<M.  The  calculation  A(M+1) 
of  the  values  for  these  functions  should  be  made  by  the  Philon  method 
[9).  Numerical  integration  of  expressions  of  type  j  a(|)co.sw£r/e 
by  this  method  is  done  for  the  number  of  points  needed  to  integrate 
only  one  function  a(c)  by  Simpson’s  method. 

The  found  sequences  of  the  functions  g  -d  are  recorded  in  the 
working  memory  of  the  machine  and  are  further  used  for  detailed  study 
of  the  antenna  characteristics.  These  characteristics  include  the 
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amplitude  and  phase  patterns  (in  any  sector  and  with  random  spacing)  , 
the  position  of  the  main  maximum,  width  of  the  pattern,  etc.  (formulas 
(7)  —  (13) ]  ,  as  well  as  the  directive  gain  and  scattering  coefficient 
(formulas  (29)  and  (30)].  Direct  integration  of  the  pattern  signi¬ 
ficantly  increases  the  time  for  determining  all  of  these  characteris¬ 
tics. 

If  we  examine  discrete  antennas  ,  the  functions  g-d  are  expressed 
in  the  form  of  sums  for  the  number  of  elements.  However,  with  a 
large  number  of  elements ,  it  is  more  convenient  to  .replace  the  discrete 
distribution  by  a  continuous  [3]. 

Random  Amplitude  Distributions 

The  technique  stated  above  is  convenient  for  computing  antennas 
by  the  method  of  statistical  tests  (Monte  Carlo  method) .  This  calcu¬ 
lation  is  reduced  to  obtaining  a  series  of  random  patterns  with 
statistical  characteristics  which  are  equal  to  the  characteristics 
of  the  original ,  and  subsequent  statistical  processing  of  these 
patterns.  In  the  field  distribution  in  the  apterture  ,  we  isolate  the 
regular  component  a^(0  (see  (1)1,  and  we  will  consider  it  to  be 
actual  (linear  phase  shift  in  the  aperture  is  not  significant) .  The 
random  distortions  of  complex  amplitude  in  the  aperture  will  be  des¬ 
cribed  by  the  functions  n  By  assigning  the  random  statistical 

nature  of  n  and  * ,  we  model  a  number  of  their  random  realizations 
(samplings)  on  a  computer  and  find  the  random  beam  patterns  which 
correspond  to  them. 

The  functions  n  and  *  are  obtained  in  the  form  of  a  discrete 
sequence  of  random  numbers.  The  sampling  spacing  is  determined  by 
the  assigned  accuracy  of  integrating  the  function  a(£). 

With  amplitude-phase  distortions,  the  vector  a(£)  is  subordinate 
to  the  generalized  Rayleigh  distribution-- n  and  x  are  independent 
normal  quantities  with  their  own  averages  ,  dispersions  and  correlation 
coefficients.  In  the  case  of  phase  errors  alone,  the  functions  n 
and  x  are  interlinked.  We  will  examine  the  latter  case  as  an 
example  ,  and  present  n  and  *  in  the  form: 

7i 13  °o  ^  (l);  *  ^  (l)  (»).  (32) 


6k 


Figure  Aa  illustrates  the  series  of  random  realizations  of  the 
function  found  on  the  computer  for  the  parameters  -^=0;  o=0.1;  p=0.1. 

Direct  study  of  the  properties  of  these  realizations  shows  a 
certain  reduction  in  their  percentage  of  randomness  (the  average  from 
several  samplings  0.080.  This  is  associated  with  the  filtering 
properties  of  the  machine  by  means  of  limitations  in  the  selected 
program. 
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Figure  5. 

The  correlation  coefficient  of  these  realizations  is  depicted 
in  fig.  Ab.  Here  the  original  curve  (33)  (solid  line)  is  plotted  for 
comparison.  Figure  5  presents  the  random  samplings  of  the  i|>  function 
with  assigned  parameters  f  =  0;  oq=0.6;  p=0.1. 

Modeling.  Random  Patterns 

In  order  to  model  the  random  patterns  ,  each  time  we  make  a  pre¬ 
liminary  computation  of  the  random  realizations  of  the  functions  n 
and  x  (in  the  examined  case  ,  only  the  functions  .  .  Then  the  obtained 
discrete  sampling  of  random  numbers  is  used  to  compute  the  functions 
gm-dm.  They,  in  turn,  are  stored  and  are  the  originals  for  analysis 
of  the  beam  pattern. 

As  an  example  ,  figure  6  presents  a  series  of  random  patterns  which 
were  computed  for  realizations  of  ij;  illustrated  in  fig.  5.  The  dotted 
line  in  this  figure  illustrates  the  avt 'age  oattern  of  the  linear 
antenna  with  corresponding  parameters. 
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Figure  7. 


Figure  7  illustrates  the  phases  of  the  random  patterns  for 
realization  of  i|i  in  fig.  5. 

When  there  is  a  fairly  large  number  of  random  realizations  ,  we 
can  construct  histograms  for  the  main  parameters  of  the  patterns. 

By  using  the  methods  of  statistical  processing  of  the  observation 
results  ,  the  statistical  characteristics  of  the  pattern  are  deter¬ 
mined.  It  is  a  priori  difficult  to  find  the  probability  properties 
of  the  patterns  in  a  whole  series  of  cases,  or  completely  impossibl 
It  is  precisely  in  these  cases  that  the  method  of  statistical  tests 
is  called  upon  to  play  a  larger  role. 


In  conclusion,  the  author  expresses  his  gratitude  to  engineer 
L.  D.  Klimenko  who  did  a  lot  oil  work  to  program  the  problem  and  the 
corresponding  calculations. 
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COMPENSATION  FOR  SYSTEMATIC  WEIGHT  DEFORMATIONS  OF  MIRROR  SYSTEM 
OF  PARABOLIC  RADIO  TELESCOPES 

B.  A.  Poperechenko  and  N.  M.  Feyzulla 

Summary 


An  examination  is  made  of  the  nature  of  the  Dhase  errors  in 
field  distribution  in  the  aperture  of  a  radio  telescope  which 
develop  under  the  influence  of  weight  deformations  of  the  main 
mirror  and  the  support  structures  of  the  irradiating  system. 

Based  on  an  analysis  of  the  normed  functions  of  individual 
components  in  the  errors  ,  a  technique  for  their  mutual  compensation 
is  suggested. 
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Introduction 

The  trend  of  developing  full-rotating  mirror-parabolic  radio 
telescopes  is  characterized  by  a  constant  rise  in  the  dimensions  of 
the  mirrors.  This  raises  the  primary  problem  of  maintaining  the  shape 
of  their  working  surface. 

At  the  early  stages  of  development  of  the  technology  of  mirror 
antennas,  the  striving  to  guarantee  accuracy  was  expressed  as  a  require¬ 
ment  that  the  deviation  from  the  ideal  profile  would  not  exceed  the 
assigned  limit  (usually  1/16  of  the  working  wavelength)  over  the  entire 
or  part  of  the  mirror  surface. 

Of  course,  this  places  excessively  strict  requirements  on  both 
the  technology  and  the  design. 

The  statistical  approach  [1,2,3]  makes  it  possible  to  consider 
more  completely  the  factors  which  influence  the  quality  indicators  of 
the  mirror  radio  telescopes. 

The  influence  of  the  statistical  factors  is  considerably  reduced 
in  the  radio  telescopes  whose  reflecting  surface  is  designed  in  the 
form  of  precisely  formed  panels  which  are  set  during  the  initial 
adjustment.  Consequently,  tne  main  factor  which  determines  the  limiting 
working  wavelength  is  the  deformation  of  the  design  under  the  influence 
of  wind  (to  a  considerable  measure  accidental)  and  weight  (systematic) 
loads. 

The  attempt  to  reduce  these  deformations  to  a  given  limit  results 
in  the  unjustified  increase  in  the  weight  of  the  metal  structures. 
Therefore  ,  in  the  practice  of  designing  large  radio  telescopes  ,  different 
methods  are  encountered  more  often  fox  compensating  for  the  deforma¬ 
tions.  These  methods  reduce  the  requirements  for  rigidity  of  the 
designs. 

There  are  known  methods  ,  for  example  ,  of  compensating  with  the 
help  of  automatic  correction  for  the  position  of  individual  elements 
in  the  working  surface  of  the  main  mirror  ,  or  counterreflector  of  the 
radio  telescope  [4, 5].  In  this  case,  a  considerable  number  of  tracking 
devices  are  used  with  fairly  cumbersome  power  actuating  mechanisms. 

The  placement  of  these  mechanisms  on  the  designs  of  the  radio  telescope 
significantly  increases  the  weight  of  its  mirror. 


It  is  therefore  more  promising  to  use  the  method  based  on  the 
mutual  compensation  for  deformations  of  different  elements  in  the 
mirror  system:  metal  structures  of  the  main  mirror  and  support 
structures  of  the  irradiating  system.  In  this  case,  the  main  mirror 
is  designed  not  for  the  minimum  of  deformations  ,  but  for  the  minimum 
deviation  of  its  surface  from  the  profile  which  is  the  most  convenient 
for  mutual  compensation  [6]. 

This  method  makes  it  possible  to  compensate  only  for  systematic 
deformations  which  develop  under  the  influence  of  the  weight  loads. 

But ,  if  one  considers  that  this  percentage  of  deformation  for  large 
radio  telescopes  is  no  loss  than  85!£,  and  with  an  increase  in  the 
dimensions  of  the  radio  telescope,  it  increases  even  more,  then  it 
becomes  understandable  why  this  method  of  compensation  is  promising 
precisely  for  large  instruments. 

This  work  examines  the  mutual  weight  compensation  for  deformations 
in  metal  structures  of  the  mirror  system  of  radio  telescopes.  It  is 
done  by  selecting  the  rigidity  of  the  support  structures  of  the 
irradiating  system  and  the  loads  on  these  structures,  and  in  certain 
limits ,  by  selecting  the  distribution  of  rigidity  between  the  design 
elements  of  the  main  mirror. 

In  this  case,  the  supports  of  the  irradf.ating  system  must  be 
designed  so  that  the  required  laws  for  movement  of  the  primary  irradi¬ 
ator  and  counterreflector  are  guaranteed  in  a  vertical  plane  in 
the  function  of  the  elevation  position  of  the  radio  telescope.  Further, 
the  indicated  longitudinal  and  transverse  shifts  in  the  primary 
irradiator  and  the  counterreflector ,  and  the  rotation  angle  of  the 
latter  will  be  called  the  compensation  parameters.  In  addition,  it  is 
suggested  that  the  elevation  correction  which  is  made  in  the  guidance 
system  of  the  radio  telesdope  be  used  as  the  compensation  parameter.  All 
the  indicated  weight  compensation  parameters  must  be  programmed. 

It  should  be  noted  that  selection  of  the  distribution  of  rigidity 
between  the  design  elements  of  the  main  mirror  makes  it  possible  to 
obtain  that  law  of  large-scale  deformations  of  its  surface  which 
promotes  the  most  complete  weight  compensation. 
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General  Correlations 

In  order  to  compute  the  compensation  parameters  ,  we  will  write 
the  equation  of  the  phase  front  [7]  in  the  aperture  of  the  radio 
telescope  with  regard  for  all  the  indicated  distortions  in  the  mirror 
system  as  applied  to  a  two-mirror  system  of  irradiation.  This  plan  is 
used  most  often  in  the  short  wavelength  region  of  the  working  range 
of  the  radio  telescope  ,  where  compensation  for  deformations  is  primarily 
needed.  We  note  that  the  weight  deformations  of  the  reflecting  surface 
of  the  counterreflector  proper  are  not  used  in  the  calculation.  They 
are  usually  insignificant. 

We  will  use  the  Cartesian  coordinate  system  with  beginning  at  the 
focus  of  an  undistorted  main  mirror  with  axis  z  which  coincides  with 
its  focal  axis  (fig.  1).  We  will  designate  the  position  of  the  phase 
center  of  the  irradiator  with  regard  for  the  weight  displacements 


Xi  -  Ja%  +  k(—  <r>  +  iVi>j) «  a  *i>  -  *<i\ 


The  weight  displacements  of  the  apex  of  a  small  mirror  will  be 
characterized  by  the  vector 

AA,=7dFy-f  A/\ 
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while  its  angular  deviation  from  the  calculated  position  under  the 

'// 

influence  of  the  weight  loads,  by  the  angle  Ac*.  f 

It  is  assumed  everywhere  that  the  weight  loads  cause  shifts  in 
the  elements  of  the  irradiating  system  only  in  a  vertical  plane.  / 

The  profile  of  the  main  mirror  with  regard  for  its  weight  defoi|na- 
tions  ( 

r  =  f{x  +  a- (*,  y ))  +  7(y  +  A.v  (*.  y ))  +  K ( <*•  ^ !  • 

where  f--focal  distance  of  the  main  mirror,  ^ 

6  ,6  ,6  --components  of  the  deviation  of  the  surface  from  parabolic 
x  y  z 

shape  under  the  influence  of  weight  loads. 

The  hyperbolic  profile  of  the  displaced  small  mirror  • ./ 


where 


*  P  +  i  [  yM  i  •  A  F,  -  A.  (/'  4-  ]  +  k  [  :M  +  AF,  + Ary M], 


2  —  __  r\  l/  }  _j_  ^  ® 

M  “  V  2  /  V  1  F(Q~F)  Tf  * 


4>--calculated  distance  between  the  foci  of  the  main  and  small  mirror, 
F--calculated  distance  of  the  apex  of  the  small  mirror  from  the  focus 
of  the  main  mirror. 

We  will  write  the  ratios  that  interlink  the  incident  and  reflected 
wave  fronts  for  both  mirrors.  We  will  express  the  single  perpendicular 
to  the  front  of  the  incident  wave  as 


and 


The  perpendicular  for  the  front  of  the  reflected  wave  is  written  as 


?  -  _ri _ y-R 

|?„-7|  and  l=T?Ts|’ 


where  Y^  and  Y-- vectors  of  the  surface  of  the  front  of  the  reflected 
wave  for  the  small  and  main  mirrors  respectively. 

The  single  perpendiculars  to  the  surface  of  the  mirrors  are 
correspondingly  equal'  to: 


yM  '  ■  *1+  (?' —  A'l)),  Y  »  #). 

The  equations  of  errors  d  in  the  form  of  wave  fronts  are 
expressed  as  a  scalar  product  of  the  vectors  for  the  deformed  and  non- 
deformed  wave  fronts : 

Here  n--perpendicular  to  the  planar  surface,  perpendicular  to  the  direc¬ 
tion  of  the  maximum  radiation  of  the  antenna.  Bearing  in  mind  the 
already  mentioned  small  elevation  correction  for  the  position  of  the 
radio  axis  of  the  antenna  A 6  ,  we  write 

The  perpendicular  to  the  undistorted  spherical  front  which  is 
reflected  from  the  small  mirror,  is  expressed  through  the  current 
coordinates  of  the  point  on  the  surface  of  this  minor: 


1  *M  "H  y.M  -T  K 


ICMi/ 


,  *.\l  +  y.M 


The  positive  values  d^  and  d  correspond  to  the  shortened  rays ,  and 
consequently,  the  leading  phase. 

The  resulting  errors  in  the  phase  front  which  take  into  consi¬ 
deration  both  the  deformation  of  the  main  mirror ,  and  the  mutual 
movements  of  all  irradiation  elements,  will  be  presented  approximately 
as  the  sum  /(-</+ </.«• 
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This  summation  surge 'cs  the  insignificance  of  the  error  in 
determining  d^  near  the  surface  oF  the  maih.  mirror  as  a  consequence  of 
ignoring  the  difference  1"  the  directions  of  the  radius  from  the 
begining  of  the  coordinates  and  the  true  beam  that  is  reflected  from 
the  small  mirror.  We  will  estimate  the  limits  of  correctness  of  this 
hypothesis . 

It  follows  from  the  equation  of  the  eikonal  for  a  uniform  medium 

that 


|  Grad /  j'->  „  1  - 1  grad, /- 1*\ 

where 

L  =  q(x,  y)  ^.mhct(v* 


grad  L--gradient  of  eikonal  in  plane  tangential  to  the  sphere  , 
s 

gradqL--the  same,  in  a  radial  direction, 
q(x  ,y) --equation  of  spherical  surface. 

By  transforming  the  right  side  of  the  equality,  we  obtain 


j  tiaA'L  |»  =  |  gradA,^  (*.  y)?~U  +  \  grad/  D(»  —  l  grad / 1)  - 

»  ( 2+  «  2  -J**- , 

\  dq  /  dq  dq 


where  we  assume 


dq 


Then  the  error  in  determining  the  distortion  of  the  phase  front 
at  the  surface  of  the  main  mirror  that  we  are  interested  in 


**)-«)—■=■- 
^  (PVK  *m  •  y,\i)  -  <?)  |  srad^i„c.r  (*■  y)  j 

j  • 

Taking  into  consideration  that  d^(x,y)  basically  has  a  large- 
scale  nature  ,  in  our  estimates  one  can  consider  with  sufficient 
accuracy  that  v' 


grad,  dM  (  v.  y)  k  grad  dM  (  v,  y). 

*  net 


For  quantitative  estimate,  we  will  adopt  the  quantity  of 
permissible  error  in  computing  the  weight  compensation  1/20  -1/30 
from  the  shortest  working  wave  of  the  TNA-1500  radio  telescope  with 
a  main  mirror  64  m  in  diameter  with  focal  distance  of  23.6  and  counter- 
reflector  5  m  in  diameter.  This  quantity  is  about  1  mm. 

Then  the  maximum  permissible  value 


!  grail </)' 


l#l-v 


10' 


This  corresponds  ,  for  example  ,  to  10  mm  deviation  in  the  phase  front 
from  the  sphere  for  1  m  of  the  displacement  on  the  sphere  surface  near 
the  small  mirror . 

A  xairly  smooth  nature  of  deformations  in  the  main  mirror  is  also 
hypothesized  in  the  further  calculations.  This  is  always  guaranteed 
in  real  designs. 

Taking  into  consideration  the  adopted  limitations  .and  omitting 
in  the  expressions  and  d  the  quantities  of  the  second  order  of 
smallness  and  the  constant  components  ,  after  all  the  necessary  calcula¬ 
tions ,  we  obtain  the  unknown  distortions  in  the  expressions  of  the 
wave  fronts; 


d  =  Ap  y  + 


x  y 

yM*.  y)—  j  fiy  (x,  »/>-*- 1*.  s*) 


H- 


4/* 


In  order  to  add  these  two  expressions  with  regard  for  the  approximation 
that  has  been  made,  we  will  express  it  through  the  same  argument  x,y. 
For  this  purpose  we  use  the  link  between  :c  ,y  and  x^,  y^  which  follows 
from  the  geometric  considerations: 

/  ©  \  -i  /  xm  '<■  y'\t  <i> 

%  _  Vm  =  l  2  ~  F )  V  1  ^  >  ~  f)  ~  '_2_  _  _fA| t 

*  y  *+*i’A  '  : 

We  obtain  a  solution  for  the  given  system  of  equations: 
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F/-(<D- 

_  **  +  </* 

Z^  +  1/1  A* 

(_©.  + 

l  4/  'J 

l  2  )  F(0~F) 

With  regard  for  the  obtained  correlations  in  the  coordinates  of  the 
main  mirror  surface 
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The  obtained  expressions  d  and  d^  may  be  used  to  compute  the 

mutual  compensation  of  the  axial  and  transverse  deformations  of  the 

main  mirror  and  the  designs  for  attaching  the  irradiating  system. 

The  design  calculation  of  the  mirror  does  not  reveal  the  three 

components  of  deformations  6x(x,y),  6y(x,y)  and  6z(x,y)  indicated 

in  the  formulas  above,  but  the  two  functions  equivalent  to  them  6  ^  ( r ,  ) 

and  <$  (r,<t>)  expressed  in  a  cylindrical  system  of  coordinates,  where 
z 

M*.  y)  =-■  br  (r.  <P)cosq>;  6y  (A-,  y)  ~  6,  (r,  <p)  sin  cp; 

M*.  ib  ="  *•  ('■.  t);  r  =  /  +  «/=';  tgq>  =  --  . 

* 

Then  in  the  cylindrical  system  of  coordinates 

—  -  -  b,  (r,  <[)  cos’  (f  —  y  br  ( r ,  <jp)  sin*  <j>  2  bz  (r,  f) 

A§rsin<jp -j - - - - - . 

•+£ 

We  transforr..  the  obtained  expression: 


i 

i 


By  isolating  in  the  expression  for  d,  the  axi symmetrical  component 
and  the  first  angular  harmonics,  we  will  obtain  the  final  expressions 
which  are  convenient  for  the  compensation  calculations: 


1  IF 

(/,  - r  A  p  siitq, 

L  ■  «  J 

Thus ,  the  resulting  error 
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H'«F 

where  AF,(p),  Ad>.(p),  AFyiP),  Aa(p),  A(l>v ( p) »  Ap(p)  --parameters  of  compensation 

and  --normed  functions  of  partial  distortions  of  the  wave  front 

for  the  corresponding  parameters  of  compensation. 

The  condition  for  complete  compensation  is  presented  in  the  form 
of  the  equality  M  to  the  constant  for  each  of  the  angles  6 .  Since 
this  condition  is  not  feasible  in  accuracy,  then  from  -$■  'we  should 
isolate  the  part  which  completely  satisfies  this  condition  with  the 
minimum  residual  part.  This  residual  in  the  final  analysis  also 
limits  the  maximum  wave  of  the  radio  telescoue. 


•'■‘i***  WMv'm 


Since  in  real  designs,  the  calculation  is  not  successful  in 
determining  and  realizing  the  compensation  parameters  with  the  necessary 
precision,  then  one  should  stipulate  the  possibility  of  regulating 
their  dependences  on  6. 

The  compensation  parameters  A/-\(p),  A<T\(p),  AO-,(p)  contain 

nonregulated  and  regulated  parts.  The  first  of  them  inevitably  arise 
as  a  consequence  of  the  nonabsolute  rigidity  of  the  design  of  the 
irradiating  system.  The  second  are  created  by  force  and  are  roughly 
regulated  through  selecting  the  rigidity  of  the  designs  in  planning  , 
as  well  as  more  accurately  by  selecting  the  weight  and  arrangement  of 
the  adjustment  loads  when  the  compensation  system  is  tuned.  The 
compensation  system  tuning  is  controlled  with  the  use  of  special 
control  and  measuring  equipment  accordng  to  the  technique  examined 
below. 


Calculation  of  the  Compensation  Parameters 

In  addition. to  the  analytical  solution  to  the  problem  of  minimizing 
the  errors  in  the  wave  front  that  requires  the  use  of  a  computer,  a 
graph-analytical  method  of  compensation  computation  is  also  possible. 

We  will  examine  it  in  the  example  of  the  TNA-1500  radio  telescope  with 
diameter  6A  m  with  focal  distance  23.6  m  (2^q=136°).  We  will  take 
two  variants  of  the  two-mirror  plan:  with  fairly  large  focal  distance 
of  the  small  hyperbolic  mirror  <5>=f=23.6  m  (distance  between  the  focuses 
and  with  fairly  small  focal  distance  5  m) .  The  latter  variant 
corresponds  to  the  distance  of  the  irradiator  from  the  apex  of  the 
small  mirror  3.5  m.  With  2^ q  =1 36 0  ,  the  parameters  of  the  mirror 
system  which  are  used  in  the  calculation  are  given  in  the  following 
table  for  both  versions. 


TABLE 


Versions 

Geometric  parameters 

of  mirror  system 

i  i 

i  1 

1 

1  F(0- 

F),  M‘ 

i 

M  | 

2 

D 

ami 

US 

O-f*.  M 

Version  I 

■ 

23,6 

39 

10 

1.8 

118  1 

21,8 

H.8 

Version  II 

5 

5,25 

1.0 

1.5 

2,5  1 

3,5 

2,5 

— 

' . . - 
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In  calculating  the  parameters  presented  in  the  table,  the 
following  correlation  was  used 


t 


I  A'/jf  (  R,\i  . 

-  V  sin  I  sill  1,  * ' cos  ^'° 
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<1.  '.a 


--  «J,  cos  M'uj  -r  4  (■•>  )  siirvj'o  j  «  Q, 


From  here  ,  with  assigned  <J> 


For  both  versions  ,  we  will  construct  graphs  for  normed  functions 
for  each  of  the  compensation  parameters  (fig.  2): 
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Compensation  for  the  weight  deformations  in  the  main  mirror  is 


achieved  by  selecting  the  scale  of  the  functions  Ai(|)  separately 
for  symmetrical  (a^  ,  A2)  and  skew- symmetric  (first  version)  (A3,An,As) 
components  as  applied  to  the  law  of  its  deformations  d(^).  As  the 
elements  which  realize  the  weight  compensation  in  the  design  of  the 
radio  telescope  ,  the  central  arrow  of  the  irradiator  and  its  fastening 
elements  on  the  main  mirror,  the  support  shape  of  the  small  mirror  ani 
the  elements  of  fastening  it  to  the  beam,  and  the  device  of  elevation 
compensation  in  the  guidance  system  are  used.  In  addition,  it  stands 


to  reason  that  measures  are  taken  in  designing  the  main  mirror  so  that 
the  functions  dg  and  dg  are  obtained  which  are  the  most  convenient 
for  compensation.  It  is  apparent  from  the  graphs  in  fig.  2  that 
depending  on  the  geometric  parameters  of  the  irradiating  system,  the 
nature  of  the  function  Ai  changes.  This  can  also  be  used  for  weight 
compensation . 

The  criterion  for  optimal  compensation  is  the  minimum  value  of 
the  standard  error  of  the  phase  front  o(B)  which  is  average  in  the 
entire  working  interval  of  the  antenna  positions  in  the  half- space: 

K 

2 

0ci>  --“Jo  (fi)  cos  p  d  p.  v 
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where 


°(P)  St  I^nucuj  (s>  P)  r  ^0U»T  (•‘>.  p)]1  ds 

--effective  (with  regard  for  the  distribution  of  amplitudes  in  the 
aperture)  standard  error  in  the  phase  front. 

The  formula  for  o(B)  designates: 

Mg)  __  aperture  distribution  of  field, 

/,m 

determined  by  expression  (3) , 

d(s,  p)  - 

doCTUT  — " 

error  in  phase  front  as  a  consequence  of  incomplete  weight 
compensation  of  the  symmetrical  and  skew- symmetrical  deformation  com¬ 
ponents  . 

The  adjustment  and  verification  of  the  system  of  weight  compen¬ 
sation  of  the  radio  telescope  uses  a  range  finder  device  which  controls 
the  displacement  of  the  datum  marks  attached  to  the  working  surface  of 
the  main  mirror.  The  two-way  device  of  the  range  finder  is  installed 
near  the  focus  or  directly  in  the  focus  of  the  corresponding  irradi¬ 
ation  circuit.  The  initial  attachment  of  the  long-range  device  is  made 
with  the  zenith  position  of  the  radio  telescope  from  the  results  of 
high-precision  geodetic  control  of  the  condition  of  its  working  surface 
and  the  position  of  the  elements  in  the  irradiating  system. 

The  functioning  of  the  compensation  system  is  checked  from  the 
change  in  distances  to  the  datum  points  with  a  change  in  the  angle  of 
the  radio  telescope  site.  In  this  case,  in  order  to  eliminate  the 
effect  of  the  variable  weight  components  of  deformations  that  develop 
under  the  influence  of  the  wind  loads,  one-sided  solar  heating,  etc. 
it  is  expedient  to  control  the  distances  to  the  datum  points  in  the 
range  of  azimuth  positions  in  limits  of  360°  for  each  fixed  angle  of 
the  site. 

For  a  one-mirror  system,  if  it  is  necessary,  compensation  of  the 
weight  deformations  can  also  be  done.  In  this  case,  the  functions 
for  compensation  for  the  axisymmetrical  deformations  Ag  and  for  the 
skew- symmetrical  deformations  Ay  are  expressed  through  displacement 
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in  the  phase  center  of  the  irradiator  Ai'z  and  Afy  by  the  formulas: 


for 

for  A/y(P) 


Distortions  in  the  phase  front  under  the  influence  of  Afz  and  Af^ 
and  deformation  of  the  mirror 

a  «  A  t,  (p)  A4  -i-  r  sin  <f  A  f,  (p)  A;  (  M  i-  d. 

The  compensation  technique  is  the  same  as  for  the  two-mirror  system. 

The  suggested  method  of  mutual  weight  compensation  of  the  main 
mirror  and  the  support  of  the  irradiator  system  is  done  by  selecting 
the  laws  of  weight  deformations  without  making  excessively  high 
requirements  for  the  absolute  rigid  metal  structures  of  the  mirror 
system.  The  calculations  confirm  the  feasibility  of  the  indicated 
requirements  and  the  obtaining  of  a  significant  weight  compensation. 

This  method  should  be  acknowledged  as  one  of  the  main  means  of 
guaranteeing  the  possible  operation  of  especially  large  radio  tele¬ 
scopes  in  the  centimeter  wave  range.  Since  the  percentage  of  weight 
deformations  as  compared  to  the  wind  deformations  rises  for  larger 
radio  telescopes  with  an  increase  in  their  diameter ,  then  the  effec¬ 
tiveness  of  this  method  for  these  radio  telescopes  must  rise.  In 
this  case,  the  optimal  movement  of  the  support  structures  of  the  irra¬ 
diating  system  can  be  guaranteed  not  only  by  selecting  their  rigidity 
and  loads  ,  but  also  by  forced  movement  of  the  elements  in  the  irradi¬ 
ating  system. 

Ir.  conclusion,  the  authors  express  their  gratitude  to  I .  V. 

Vavilova  for  detailed  and  useful  discussion  of  the  article  manuscript. 
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INPUT  IMPEDANCE  OF  SYSTEM  MADE  OF 


IN TEES EC TINS 


LOW-SET  VIBRATORS 


V.  P.  Kul’tsep 


Summary 


The  Hansen-Eekerly  method  is  vised  to  determine  the  input  impedance 
of  a  system  made  of  symetrical  vibrators  with  sinusoidal  current 
distribution  which  intersect  at  the  power  point  and  are  suspended 
at  a  low  altitude  above  the  earth. 

A  single  low-set  vibrator  is  examined  as  a  particular  case. 

Results  are  presented  from  measurements  of  the  input  imoedance  of 
the  system  of  intersecting  vibrators  in  the  short-wave  and  medium- 
wave  ranges . 


Introduction 

In  determining  the  complete  input  impedance  of  a  low-set  symmetrical 

vibrator,  the  latter  is  generally  represented  by  equivalent  segments  of 

a  long  line.  The  linear  parameters  of  the  line  are  found  by  analyzing 

a  line  of  infinite  length  with  the  same  geometry  of  the  cross  section. 

This  calculation  provides  values  of  the  reactive  part  of  the  complete 

input  impedance  which  are  satisfactory  for  practice.  The  results  of 

this  calculation  show  that  the  input  impedance  of  a  low-set  vibrator  is 

mainly  determined  by  the  linear  successive  impedance  R-.  which  lies  in 

0  ^  o 

limits  from  30it  /\  Ohra/mf°r ground  with  infinite  conductance  to  60n  /X 

Ohm/m  for  purely  dielectric  .croundf  1  ,2 ]  . 

The  generalized  method  of  How  makes  it  possible  to  obtain  a  constant 
component  of  the  input  impedance  which  does  not  depend  on  the  frequency 
[1].  Dipole  [3]  and  multipole  components  are  also  known  which  are 
proportional  to  the  square  and  higher  degrees  of  frequency. 

All  the  components  can  be  taken  into  account  by  introducing  a  cer¬ 
tain  equivalent  linear  impedance  R^  which  for  short  conductors  may 
significantly  exceed  [4] ,  or  by  computing  the  input  impedance  by 
one  of  the  strict  methods.  Two  methods  are  quite  strict  on  the  con¬ 
dition  that  the  current  distribution  is  assigned:  the  directed  emf 
and  the  Hansen- Bekerly .  These  methods  are  equivalent  in  the  strictness 
of  general  expressions.  However,  in  computing  the  near-earth  antennas 
by  the  first  method,  the  emf  which  is  directed  by  a  field  reflected 
from  the  earth,  is  only  taken  into  account  approximately,  especially 
for  the  induction  components  of  the  field.  In  particular,  the  con¬ 
stant  component  of  losses  in  computing  by  the  directed  emf  method 
generally  drops  out  of  the  solution.  At  the  same  time,  this  component 
has  a  significant  value  for  the  electrically  short  antennas. 

The  purpose  of  this  work  is  to  determine  the  input  impedance  of 
a  system  made  of  N  symmetrical  vibrators  2L  long  and  d  in  diameter 
which  intersect  at  the  power  point  at  angle  and  which  are  suspended 
at  altitude  h  above  the  surface  of  homogeneous  ground  with  permeability 
e  and  conductance  a  using  the  Hansen-Bekerly  method.  The  input  impe¬ 
dance  of  a  single  vibrator  is  found  in  a  particular  case. 
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The  llnnson-Bekerly  method  has  already  been  used  to  determine  the 
input  impedance  of  a  single  vibrator  ,  however  ,  the  analysis  either 
did  not  take  into  consideration  the  most  significant  induction  compo¬ 
nent  of  losses  [5]  or  was  limited  to  the  case  of  a  low-conductiig 
ground  16)  , 


Statement  of  the  Task 

We  will  examine  an  idealized  system  made  of  2N  flat  current 
sectors  with  radius  L  located  at  altitude  h  (fig.  1).  We  assume  that 
the  average  width  of  the  sector  equals  the  diameter  of  the  vibrator 
( d =<}> qL )  .  We  make  an  axis  of  the  cylindrical  coordinate  system  p  ,<{>  and 
z  through  the  intersection  point  of  the  sectors. 


It  is  known  [7]  that  the  stream  of  active  power  which  passes 
through  the  plan  located  above  the  emitter  with  assigned  current 
distribution  ,  precisely  equals 

=  60a* k  V  V 

s— 2  0 


0) 


The  stream  of  active  power  that  passes  through  the  plane  which 
is  located  below  the  emitter,  but  above  the  earth's  surface: 


In  these  expressions 


p  d  (>  d  (J  dr 


(3) 
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--coefficients  (jjL  expansion  of  the  vector  potential  of  the  assigned 
distribution  Ksifj  the  current  density  vector  i  in  the  system  of  ortho¬ 
gonal  vector  functions  A±  which  form  a  discrete  spectrum  of  azimuth 
harmonics  with  different  n  and  noncontinuous  spectrum  of  waves  with 
different/ angles  2V=arc  cos  1/k  (including  the  complex  with  1  >  k)  be¬ 
tween  thet  Wave  perpendicular  and  the  plane  of  the  horizon.  The  re- 
maining  ck/i/l ^nations  have  the  following  meaning: 


**  ~  m\  t  m  — 

*r  ~  % 

Hi  .  .  Hi.  ’  1  r  M. 


IW  ;  Wj 


t  m  -j-  m,  * 


m  m, r«  £  -i  ip. 

P  -  COo/.,  2.-1/)., 

X--wavelength. 

We  assume  that  the  vector  of  current  density  has  only  a  radial 
projection  which  is  assigned  by  the  product  of  two  functions 

if  ™  Vi  ((*)/*  (T).  (5> 

where  lQ--ef f ective  value  of  the  current  at  the  inlet  of  the  emitting 
system. 

Since  i^=iz=0,  then  we  need  the  projections  Xgnonly  on  the  p  axis: 

/ 

(6* 

k  d  i> 

Then  the  coefficients  fgn  can  be  written  as  follows: 
fe  — <«•*£*•  W*  <*> 


where 


<D„  -  f  /s,w)v 
6 
L 

Tsn  Jr.il 


7';,,  /.O'!  !>*!’  «!'• 

a  i» 
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We  further  assume  that  the  normed  distribution  of  -nt  density 

along  the  vibrators  has  the  appearance  of  a  standing  wave : 


/.( P)- 


1  sin  (L  —  p) 
p  sin  kxL 


(P\ 


where  k^=nk,  n-- shortening  of  the  vibrator. 

This  case  does  not  take  into  consideration  the  non-inphase  nature 
and  damping  of  the  amplitude  of  current  distribution  along  the  vibra¬ 
tor  which  are  associated  with  radiation.  However,  these  effects  are 
comparatively  small  for  electrically  short  vibrators. 

Normed  azimuth  distribution  of  the  current  density  for  an  odd 
number  of  vibrators  is  assigned  in  the  form 


where 


(+  l/.V«fowith,h'<<r<'h2- 
h  (<f)  -  —  l/A'  <r,with*  +  <p!i  <  <p  <  n  +  ‘P'** 

vO  outside  current  vectors 

<Fu  *= » <Pi  “  <Po/2,  fin  =  » 


(14) 


i- -number  of  vibrator. 

In  computing  (10)  with  regard  for  (14)  we  have  for  odd  azimuth 
harmonics 


sin  Nn  tpi/2 
sin  n  <p,/2 


(15) 


TJ’or  n  even  $  =0. 

n 

It  can  be  shown  that  (15)  is  correct  even  for  an  even  number  of 
vibrators . 

In  order  to  compute  T2n  we  use  the  presentation  of  Besselian 
functions  through  the  producer 

e(v-im„7  =  vj4(z)r*.  (16) 

Assuming  the  argument  of  the  producer  function  W=se^x  ,  we  have  the 
known  correlation 

M 

sin(?sin  \|)  —  2  ^^((J)sin(2i  +  1)  ij-,  (17) 

(-0 
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while  assuming  W=iex ; 


?in  (:rh /)---  2  V(-l)03|fl(?)ch(2/+  I)*. 

tu*0 


By  transforming  the  numerator  of  (13)  with  l>k^  with  the  help  of 
(17)  ,  and  with  l>k^  with  the  help  of  (18)  ,  and  by  using  the  known 
formula 

f  -7  •'„«(«>.  U9) 

0 

we  obtain  the  following  series  for  T2n-* 

Tin  =  “  Yi  ^<+l  <20^ 


where  -j/+i  = sln  (2<  +  1 )  »J>,  \j>  =  arc  sin  kjl,  if  l>k  and  t2i+1=  (— l)‘ch(2/+ l)x,  . 

’,“,"[T+v/Tt")’-1]  , « x<kr 

We  similarily  obtain 

T“  “  2  tT  E  ‘ 1 '  («•>  +  2  jj  “■/  X 

L  t-o  /"i 

X?!'-»'W, +.<“•>]•  (!l) 

i-0  -* 

where  U2j-cos2ji}>  with  l>k^  and  U2j=(-l)Jch  2j  with  l<kj . 

The  substitution  of  (20)  and  (21)  into  (8)  and  (9)  and  further 
in  ( 1 )  —  < 2 )  provides  an  accurate  solution  to  the  set  task  in  the  form 
of  the  sum  of  single  integrals  which  can  be  numerically  found. 

We  now  obtain  the  calculated  ratios  after  making  certain  permissible 
approximations . 


In  computing  R^.  ,  the  integration  variable  1  does  not  exceed  k. 

If  the  electric  length  of  the  vibrator  k-^1  is  not  more  than  115°  ,  then 
the  argument  of  the  Besselian  functions  in  (20)  and  (21)  does  not 
exceed  2.  Following  [5,6]  we  replace  Jn(lL)  with  the  first  term  in  the 
exponential  series  expansion: 
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In  addition,  since  in  the  region  which  is  important  for  the  integra¬ 
tion,  k^/l=ch  *>>1,  then  the  approximation  is  permissible 


By  using  (22)  and  (23)  ,  we  have 

I/21 1  =* ,  fs,  -1 1  <7i  •  I  hi  \  —  !  f 3,  -1  i  = 


where 


_*,L  .  [.  __  (kiL)±  Jm  (*1  L)* 
2  sin  [l  12  1  360 


and,  taking  into  consideration  only  the  terms  with  |n|-l,  we  obtain 
for  the  radiation  impedance 

; :  10<74£2/.2 1 1  -|-  -3~  f  I  |  i2  |J  2cos  2 nth  Re  aa  —  2si«  2 mh  Im  a*  4-  • 
l  *o  L 

4-  (1 3tj  Is  —  2cos  2 mh  Re  a3  f  2sin  mh  Im  a3)j  d  jj .  (25) 

We  will  also  present  the  reflection  coefficients  in  the  form  of 
a  series 

at«-l  +  -l7T- ^(T)+¥(T)-'5rlTj  t  (26) 

-mrl"  £<*-)-  -];  « 

where  r^p+\q*=?\  /  Y *—  1. 

In  computing  the  radiation  impedance ,  one  can  be  limited  to  two  terms 
in  (26)  and  one  in  (27).  Then  computation  (25)  yields 


Rv  -  lO?! k'L 


■w.-  i-l- + ( _i±s4i_y + in  -  ^ + 

|  U  - 1|  \  |  r|  -i- 1  +  2Rc  V  1  '  *  *■  5 

+  Uellg(,L_2^_|i  +  2ii)]j.  (2 
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where  l^Zkh. 


If  h'l^l,  then  with  height  of  the  suspension  approaching  2ero , 
(28)  becomes 

(29) 

i.e.  ,  for  short  and  thin  (q^-*8)  vibrators,  formula  (29) 

has  the  appearance  of  Rudenberg*s  formula  with  a  numerical  coeffi¬ 
cient  that  is  reduced  because  of  the  influence  of  the  ground. 

In  another  extreme  case,  with  h-*«  and  T>>1 

(30> 

In  particular,  for  a  thin  half-wave  vibrator  R  =73.5  Ohm. 

u 

Consideration  for  the  subsequent  terms  in  the  series  (22)  ,  (23)  , 
(26)  and  (27)  ,  as  well  as  the  terms  with  jn|  >1  yields  multipole 

4 

components  of  radiation  resistance  which  are  proportional  to  (kL)  , 
(kL^)...,  whose  quantity  is  small  for  short  vibrators. 

Loss  Impedance 

•The  dipole  component  of  losses  for  radiation  into  the  ground  of 
waves  which  have  in  the  ground  an  angle  of  inclination  of  the  wave 
perpendicular  in  limits  from  ^=90°  to  lS^=arc  cos  is  determined 

by  the  terms  with  |nj=l  in  the  first  term  of  (2).  By  computing  under 
the  same  conditions  ,  as  in  the  computation  of  R^  ,  we  have 


The  results  of  the  computation  according  to  formula  (31)  coincide 
with  the  results  obtained  in  [5]  by  numerical  integration. 

For  ground  with  ,  formula  (31)  becomes 

R,  -21  Jq***!*:?-'*.  (32) 

The  losses  associated  with  the  dipole  waves  are  emitted  into  the 
ground  with  slipping  angles  in  the  ground  from  1^=)^  to  ^=0.  They 
are  determined  by  the  terms  with  { n | =1  in  the  second  term  (2)  in  the 


interval  of  integration  k<l<k/e  : 

k  !*• 

#//  -  210.1-  /;  j  l(  Im  >  ,  U]  ,*  l'11  >»)  **_!."  •  M 

k 

Since  in  this  case  l>k,  then  a  larger  number  of  terms  is  needed 
in  the  series  (26)  to  approximate  c^.  In  the  interval  0<v< /e-1  four 
terms  provide  sufficient  accuracy: 


Im  a,  -  2p  V  +  2 pq  »*  —  (p3  +  3 pq!)  >\  (34) 

where  '"I  /*— *2/*- 

For  Imag  in  the  interval  0<v< |7-1 1 ,  the  following  expression 
yields  satisfactory  accuracy: 

Im  a3  s-  2*  Re  (s  /  \r  —  1 ) ,  (35) 

and  in  the  interval  ,/S<v<l—l- 

Im  a,  =  0,8  /  v  j  |/7—  1  j.  <36> 


To  simplify  the  computations,  we  assume  that  (36)  is  correct  for 
the  entire  interval  0-1 /e-1 ,  which  yields  a  result  which  is  exaggerated 
by  several  percentages.  Computation  of  (33)  with  regard  for  (34)  and 
(36)  ,  bearing  in  mind  that 


yields 


!  R„  -  15^4*is  {(2 P  +  0,s;r —  1  r',?)|  I  -  (1  +  Oe  W  4- 
i  +  2pql2~(2+2t  +  (l)<T,]V-i?-3pq,)X 

X  (6— {6  +  6/  +  3P)  (38)  ^ 

where  /=£)/  e— 1. 

In  contrast  to  the  expression  obtained  in  [6] ,  expression  (38) 
yields  a  final  result  even  with  h-*P.  The  difference  is  explained  by 
the  fact  that  the  approximation  adopted  in  [6]  Ima^  with  l=k  yields  a 


divergent  result ,  while  in  reality  Imn^O  with  l->k. 

For  a  short  vibrator  located  on  the  surface  of  the  ground  without 
losses,  on  the  condition  that  e>>l,  formula  (38)  becomes 

=  17,3ASZ.5 1  T-T  (30) 

A  similar  formula  was  obtained  in  [3]  for  the  input  impedance  of 
an  infinitely  small  dipole  lying  on  the  interface  of  the  dielectric 
and  the  air.  The  numerical  coefficient  in  [3]  is  larger  than  in  (39), 
as  a  result  of  the  fact  that  the  current  distribution  is  assumed  to 
be  rectangular,  as  here,  and  not  sinusoidal. 

For  low-conducting  ground,  the  component  RII  mainly  determines  the 
input  impedance.  However,  in  the  case  of  a  semi-conducting  ground, 
the  induction  losses  are  dominant.  They  are  determined  by  the  remaining 
part  of  the  second  term  in  (2). 

We  will  designate 

=  240V ?|/Si’  lm >.  .  (40) 

ft,,-  240.-.’*  (41) 


As  is  apparent  from  (20)  and  (21)  ,  the  value  of  any  of  the  functions 
Tgn  is  negligible  until  the  argument  1L  exceeds  the  least  of  the 
indices  of  the  Eesselian  functions  which  are  included  in  the  formulas. 

In  the  transition  of  the  argument  through  the  value  lL=n,  the  value 
of  the  function  drastically  rises,  and  with  a  further  rise  in  the 
argument,  it  fluctuates  around  a  certain  average  value. 

By  using  the  threshold  nature  of  change  in  the  functions  Tgn ,  we 
replace  them  by  exponential  functions  equal  to  zero  at  lL<n ,  and  with 
1L>  n ,  equal  to  the  average  value  of  Tgn  during  the  period  of  change 
in  In(lL): 


0*  T  ~  .1*  1 

IL  ‘  Tan  ~  4L  r  <«  ’ 


3 

I 


1 

j 

3 

J 


T** 

where  tis^kxL  c'.g  kxL. 


(42) 
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The  rapidly  changing  function  |<>n|  is  also  nporoximatcly  replaced 
by  its  average  value,  equal  to  A  in  the  interval  0<n< n  /2N<>^,  equal  to 
A/N  in  the  interval  n^n*^*^/^  and  close  to  0  with  n>n^. 

For  Imag  the  relative  error  on  the  order  of  57,  is  yielded  for  the 
semi- conducting  ground  by  the  approximation: 


Ini >, « c-M5./V'ni  —  F. 


(43) 


For  the  approximation  Imaj  we  use  the  asymptotic  value; 


Im 


_ _ 
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Since  the  terms  in  the  sum  in  (AO)-(Al)  differ  little  from  each 
other,  we  will  replace  the  summation  by  an  integration  according  to 
the  formula 

y 

/<*>»  ~  f  m**>  <45> 

•-Hi  2 


where  II — summation  spacing  equal  in  our  case  to  two.  We  will  further 
change  the  order  of  integration 


Ru,  =  240* 


*,/i  It  •t/L  IL 

~~  f  l)dn  +  y  f  dl\f(n,  l)dn 

l/L  «  l/N  *0 


,  (46) 


where 


.f 


a  «  2A  +  0.226s.  &  «  2li  +  0,495?.,,  -  /.  /  ]  ’FT. 


In  the  expression  /l+v*  ,  we  ignore  the  1  which  yields  a  slight 
error  for  short  vibrators  (kL<l)  since  the  lower  limit  of  integration 
for  1  is  considerably  greater  than  k.  For  longer  vibrators,  this 
ignoring  yields  a  somewhat  exaggerated  value  of  the  losses.  Finally, 


»  - (*>-  [Ei  (-  -m ■)  ■ - Ei -&)  - 

-BH)+E,(-i)]+v[H(-*;-)- 
-b(-*a)-o(--j-)+h(- ;-)]>. 


(47) 
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where  ¥*“ !*  </?+'/! A*. 

The  component  R1V  is  found  in  exactly  the  same  way: 


Kn.  - 


kL  (t 


♦  WW"  TSS-)-BK)lf 

J- 1  |Ei(-r) 


where  <74  «=  <7*  +  (A’4/-)J  </j  f-  (#i£)4/3. 


Discusssion  of  Results 

The  component  Rjjj  is  negligble  for  a  single  vibrator  shorter  than 
A4ahd  also  if  the  height  of  its  suspension  significantly  exceeds  A^. 

If  the  length  of  the  vibrator  is  much  greater  than  the  wavelength  in 
the  ground,  then,  assuming  that  L>>b,  and  by  using  the  approximation 
Ei(-x)=ln  1.78x,  we  find 


«  592—  <73  jl  n  l,275[Ei  -  Ei 


For  a  thin  vibrator  .vrtiose  height  of  the  suspension  and  the  diameter 
(width)  are  much  smaller  than  Aj  ,  formula  (50)  can  be  simplified  even 
more: 

R,n  =  2lVs~  •  (50) 

The  last  cofactor  in  (50)  is  the  linear  impedance  of  successive 
losses  of  the  conductor  of  infinite  length.  Consequently,  the  compo¬ 
nent  Rjjj  is  mainly  determined  by  the  successive  losses  of  the  vibrator. 

The  cofactor  q  depends  on  the  shape  of  the  current  distribution. 
For  a  short  vibrator  with  triangular  distribution,  q3  “1/3.  For  a 
half-wave  vibrator,  q^  *1. 

The  component  for  a  single  vibrator  with  B>>e  and  h+0: 

(5I> 

and  is  known  from  the  theory  of  grounding  electrodes  as  the  impedance 
of  spreading  for  a  dipole  made  of  uninsulated  conductors  placed  on  the 
surface  of  the  ground.  With  £>>e  and  k^L,  <<]  ,  this  component  does 


-  ..gyjgifegjfSaB 


not  depend  on  the  frequency,  and  consequently,  is  a  constant  component 
of  losses. 

Thus ,  all  the  components  of  the  input  impedance  of  the  vibrator  are 
found  by  a  single  method.  In  this  case,  in  contrast  to  those  known  in 
the  literature,  dependences  are  obtained  for  all  the  components  on  the 
height  of  the  suspension,  the  width  of  the  vibrator  and  its  length  which 
are  also  correct  for  the  limiting  values  of  these  parameters. 

The  components  and  R^  for  low- set  antennas  are  generally  smaller 
than  the  others. 

The  component  Rjj  prevails  when  the  length  of  the  vibrator  is 
smaller  than  the  length  of  the  wave  in  the  ground.  If  the  length  of 
the  vibrator  is  considerably  greater  than  then  the  component  RIII 
becomes  dominant.  The  component  R^y  has  a  noticeable  value  where  the 
height  of  the  suspension  is  commensurate  with  the  width  of  the  dipole 
even  above  dry  ground  when  e  and  8  have  one  order.  Generally,  with 
an  increase  in  the  height  of  the  suspension,  the  constant  component 
diminishes  considerably  faster  than  the  others. 

The  components  Rj,  ,  Rj  and  Rjj-  depend  on  the  width  of  the  vibrator 
as  *1*1  {2L;d)~  sin*  J/24,  i.e.  ,  very  weakly.  The  dependence  of  Rjjj 

on  the  vibrator  width  has  a  completely  different  nature. 

Calculations  according  to  formula  (49)  show  that  the  input  impedance 
of  a  fairly  long  vibrator  drops  as  its  width  increases  as  if  it 
consisted  of  N  mutually  isolated  vibrators.  The  number  N  is  determined 
by  the  number  of  half-waves  in  the  ground  which  is  laid  along  the 
width  in  the  middle  of  the  vibrator  arm. 

The  fact  that  parallel,  low-set  vibrators,  separated  by  a  half¬ 
wave  in  the  ground,  can  be  considered  isolated,  is  known.  Such  an 
assertion  is  proved  for  vibrators  which  intersect  at  one  point  with 
the  help  of  formula  (47)  with  N#l. 

In  fact,  for  vibrators  with  length  much  greater  than  X3  ,  the 
component  Rjjj  dominates.  Then  the  dependence  of  the  input  impedance 
on  the  distance  between  the  ends  of  the  dipoles  ,  computed  according  to 
formula  (47)  with  h<<.\^  has  the  appearance  shown  in  fig.  2. 

It  is  apparent  in  fig.  2  that  the  vibrators  can  be  considered 
completely  isolated,  then  D/x^  >3.  In  this  case,  the  input  impedance 
of  the  system  R^x  ^  is  N  times  smaller  than  the  impedance  of  a  single 
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vibrator  K_„  , 
&x  1 


If  the  distance  between  the  ends  of  the  vibrators  equals  the  length 
of  the  wave  in  the  ground ,  then  the  input  impedance  of  the  system 
roughly  equals 


R 


uM 


/ 1  ;  Ai‘v“  1 

aTY  1  0,4  "it' 


(52) 


The  presence  of  a  second  term  in  the  parentheses  is  explained  by 
the  coupling  impedance.  For  the  case  of  two  intersecting  vibrators, 
whose  ends  are  separated  by  ,  the  coupling  impedance  increases  the 
input  impedance  by  20Z  as  compared  to  the  case  of  complete  isolation. 

With  an  increase  in  the  number  of  vibrators  ,  the  coupling  impedance 
increases  somewhat ,  and  as  a  result ,  the  input  impedance  of  the  array 
made  of  a  large  number  of  vibrators  whose  ends  are  separated  by 
is  roughly  equal  to  R6*1^N* 

The  condition  can  be  considered  a  criterion  for  designing  the 

array  made  of  intersecting  low- set  vibrators  above  semiconducting  ground. 

Examples 

We  will  take  as  the  first  example  the  dependence  already  examined 
in  tAJ  and  [6]  of  the  input  impedance  of  a  single  half-wave  vibrator 
on  the  height  of  the  suspension  with  A =5  m. 

In  figure  3,  curves  1  and  2  show  the  calaculated  dependences  which 


were  taken  from  (4).  Curve  1  is  calculated  by  the  method  of  directed 

_  o 

emf  for  ground  with  c**10,  op10  mho/m.  Curve  2  was  computed  by  the 

•  9 

Hansen- Bekerly  method  [6J  for  ground  with  c*7,  o*2.8  x  10  mho/m,  but 
without  consideration  for  the  Rjjj  and  components.  Curve  3  was 


Figure  3. 


Figure  A. 

computed  with  the  help  of  formula  (28),  (31),  (38),  (A7)  and  (A8)  for 
ground  with  e=10  and  o=10”^  mho/m.  The  shortening  was  determined  by 
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Figure  5. 


the  technique  presented  in  [2].  The  dotted  line  shows  the  calculated 
curve  obtained  on  the  assumption  that  n*l. 

Comparison  shows  that  with  h>0,  IX  ,  the  course  of  all  the  curves 
is  roughly  the  same,  and  curves  1  and  3  almost  coincide.  However, 
with  low  heights  of  the  suspension,  the  curves  significantly  differ, 
and  calculation  hy  the  formulas  presented  here  yields  finite  values  of 
the  input  impedance.  This  really  occurs  in  practice. 

We  will  examine  as  a  second  example  the  dependence  on  the  frequency 
of  the  input  impedance  of  a  three-beam  .  antenna  (N*3)  located  on 
the  surface  of  the  ground  with  e«5,  o=10”^.  The  calculation  was  made 
according  to  the  five  formulas  indicated  above  for  the  case  h«0.  The 
length  of  the  vibrator  arm  L=4.2  m,  diameter  of  the  conductor  d»10  tnm. 
Shortening  was  assumed  to  be  equal  to  n=| /T  | .  The  aperture  angle 
between  the  extreme  beams  is  2^=60°, 

In  figure  4 ,  the  solid  line  shows  the  calculated  dependence  of 
the  input  impedance  on  the  frequency  in  the  range  2-12  MHz.  The 
points  on  fig.  4  show  the  values  of  the  input  impedance  measured  in 
the  summer.  The  prototype  was  made  from  RK-3  cable  placed  on  the 
surface  of  the  ground  (dry  sand)  and  slightly  sprinkled  with  sand 
from  above. 

As  the  last  example  we  will  examine  the  dependence  of  the  input 
impedance  of  a  two-beam  antenna  on  the  angle  between  the  vibrators  <t>^. 

The  calculation  was  made  for  the  system  made  of  two  vibrators  2L*» 

2  x  150  m  long  placed  at  a  height  of  10  mm  above  the  ground  with 
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frequency  of  360  kHz. 

Shortening  according  to  the  experimental  data  was  assumed  to  be 

equal  to  n=1.37.  The  calculation  results  are  shown  in  fig.  5.  Curve 

-3 

1  was  obtained  for  ground  e=7,  o=10  moh/m,  curve  2  for  ground  e=10, 

_  *3 

o  =1 0  moh  /m . 

The  points  show  the  measured  values.  The  vibrators  were  made  of 

2 

BPVL  conductor  2.5  mm  ,  freely  placed  on  the  surface  of  the  ground. 

The  measurements  were  made  in  winter. 

We  note  that  the  dependence  can  be  used  to  measure  the 

conductance  of  the  ground  in  the  medium  wavelength  region. 

Conclusions 

The  Hansen-Bekerly  method  makes  it  possible  to  obtain  all  components 
of  the  input  impedance  of  low- set  vibrators  without  exception.  In 
contrast  to  ethers,  this  method  yields  a  final  result  which  agrees 
with  the  experiment  even  for  heights  of  the  suspension  close  to  zero. 

Mutual  impedance  between  the  intersecting  low- set  dipoles  is  low 
if  the  distance  between  their  ends  exceeds  the  wavelength  in  the  ground. 
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